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ABSTRACT 


This report extends previous investigations of 24-hour satellites by 
considering the long-term drift of such satellites with orbits having any 
inclination but small eccentricity in an earth gravity field through fourth 
order. 

It is demonstrated that the drift regime for the inclined- orbit 
satellite is the same as for the equatorial satellite, modified only by an 
"inclination factor" distinct for each relevant gravity harmonic and in 
each case less than unity for non- zero inclinations. 

In particular, it is shown that for the circular orbit satellite 
through fourth- order earth gravity, only the longitude gravity harmonics 
J 22 , J 3i > J 33? J 4 2 > and J 44 can have long-term (secular) effects on the 
orbit. Except for very high inclinations, the 24-hour drift regime, as 
influenced by presently determined earth gravity, is dominated by the 
effect of the J 22 harmonic (that associated with the ellipticity of the 
earth’s equator). 

For the longitude gravity harmonics, J 31 and J 42 , certain "non- 
resonant" inclinations are shown to exist. Twenty-four-hour satellites 
having orbits with these inclinations will experience negligible long- 
term accelerated drift due to either J 31 or J 42 . 
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THE DRIFT OF AN INCLINED-ORBIT 24-HOUR SATELLITE 
IN AN EARTH GRAVITY FIELD THROUGH FOURTH ORDER 


by 

C. A. Wagner 

Goddard Space Flight Center 


INTRODUCTION 

In recent years there have been many studies of the perturbations of a 24-hour satellite due 
to the high- order gravity field of the earth as well as the gravity fields of the sun and moon 
(References 1, 2, 3, 4, 5 and 6). Since the period of this satellite is almost the same as the earth’s 
rotation period, its mean geographic position may be virtually unchanged for many months. It is 
evident that even very small longitude mass inhomogeneities within the earth may be able to cause 
appreciable longitude drift of this satellite, which can respond in the same direction to these dis- 
tant mass concentrations over a long period of time. This phenomenon is often referred to as 
’’resonance” in the literature. While use of this description can cause confusion if one visualizes 
the action only in geographic space, it is an accurate one for the action in inertial space and it will 
be retained. 

A virtually fixed station in the sky with almost every place on earth permanently visible from 
one of three such equatorial stations is an obvious advantage to worldwide communications and 
navigation. Previous studies of the equatorial 24-hour satellite (References 2 and 5) have shown 
that the long-term east-west station- keeping requirement for this satellite is controlled by second- 
order earth- longitude gravity which is associated with earth- equatorial ellipticity. 

The first operational 24-hour communications satellite, Syncom II, launched in July 1963, has 
an orbit with an inclination of about 33°. The study in Reference 6 shows the drift regime of such 
a satellite in a second-order longitude- gravity field to be essentially unchanged from that of the 
equatorial satellite in the second-order field. Syncom II has now (Spring, 1965) drifted almost 
around the earth in over a year. Comparison of observations on its drift in this time with the 
gravity drift theory of Reference 6 has already revealed the ellipticity of the earth’s equator to a 
greater precision than ever before (References 6 and 7). The advantages of the use of freely drift- 
ing 24-hour satellites for geodetic purposes have thus already been dramatically demonstrated. 

This study is designed to provide a simple theoretical framework for the use of such satel- 
lites of any inclination for basic investigations of the complex gravity field of the earth to high 
order. 
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SUMMARY OF PREVIOUS INVESTIGATIONS 

In Reference 5 it was shown that in the vicinity of a momentarily stationary longitude K 0 the 
longitude-drift acceleration of a 24-hour, circular orbit, equatorial satellite in an earth gravity 
field to fourth order is given by 

X o = Y 22 sin2 K -K aa ) + F 31 sin(k 0 -k 31 ) + F 33 sin 3(\ 0 - k 33 ) 

e s ^ 


where 


+ F 42 sin 2(\ 0 -\ 42 ) 


F 44 sin 4 


(^0 


rad/(sid. day) 2 


(1) 


F 2 2 “ a 2 ( R o/ a s) 2 J 22 * 

s 

F 3i * 2a 2 ( R o/^s) 3 J3i » 

45 fJL z 

F 33 " “ 2 ( R o/ a s) 3 J 3 3 • 

~15^ e 

F 4 2 “ a 2 ( R o/ a s) 4 J4 2 * 

s 

42°m e 

F 44 ' a 2 ( R 0 /^s) 4 J44 

s 


(2) 


The constants are defined as follows: aj e is the rotation rate of the earth (2 ? t radians/sidereal 
day); a s is the "synchronous" semi- major axis of the satellite; is the Gaussian constant of the 
earth; R 0 is the mean equatorial radius of the earth; J nm and \ nin are earth gravity constants in the 
spherical harmonic expansion of gravity for a general mass distribution. The J nm (m / 0) for the 
earth are very small dimensionless constants, of the order of 10~ 6 or less. The j nm give the 
magnitude of longitude- dependent gravity forces. The \ nm give the phase angle of the dependence. 
For a truly spherical earth, every J nin would be zero (Reference 8). The F nm sinm(A. 0 “-\ nm ) terms 
in Equation 1 are the earth- gravity longitude perturbation forces (on the equator) for the particular 
harmonic, nm. Equation 1 is a statement of the proportionality of these longitude driving forces and 
the resultant drift acceleration for the 24-hour equatorial satellite. It is not a self-evident statement 
because the coordinate system r, <£ is not inertial but is rotating at the earth rate with respect 
to inertial space. 
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In Reference 6 the orbit- averaged longitude drift of a 24-hour, inclined circular- orbit satel- 
lite in a second-order longitude- dependent earth field was found to satisfy the equation 

k = A 22 sin 2 ( K -^22 ) * (3) 


where 


A 22 = -72t t 2 (R 0 / a s) 2 J22 ( C ° S 2 + *) rad/ ( sid * da y) 2 . 

and the gravity constant k 22 locates the longitude of the major axis of the hypothesized elliptical 
equator of the earth. Equations 1 and 3 are compatible since A 22 is equal to 

3(2*) 2 F 22 


for i = 0, noting that 


defines the semi- major axis of the circular- orbit earth satellite with a period of exactly one day. 
The strict derivation of Equation 3 applied only to a 24-hour satellite in the vicinity of a stationary 
configuration whose ascending node (or mean daily geographic longitude) was located at k. How- 
ever, even for the long period and wide longitude excursion of libratory drift which follows from 
the regime established by the '’pendulum equation,” Equation 3, the near- stationary conditions 
assumed in its derivation are not seriously violated (Reference 6). Thus, for the entire libratory 
regime, Equation 3 can be expected to hold sufficiently well with\ representing the daily geographic 
position of the ascending node. In fact, numerically integrated Syncom II trajectories in the 
presence of sun, moon and earth zonal gravity have shown that the drift regime established by the 
simple formula, Equation 3, with k interpreted as the daily position of the ascending node of the 
satellite, is a fully adequate representation of the actual long-term drift within the small ’’noise 
limits” of these higher order gravity perturbations (References 6 and 7). These trajectories have 
included some with eccentricities as high as 0.0012. Drift rates have been as high as 0.8 degrees/ 
day, which is twice the maximum allowable for this 24-hour satellite in a libratory regime (Refer- 
ence 6). R. R. Allan* presumes that Equation 3 may fairly represent the 24-hour drift regime 
for eccentricities as high as 0.3 if k is taken as the mean longitude of ascending and decending 
nodes which, for perigees away from the equator, is a better measure of the mean daily geographic 
longitude than either nodal position alone. 


•Private communication. 
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At any rate, for a sufficiently wide specification of orbital elements to be practical, Equation 3 
has been shown to establish the correct second-order longitude drift regime with K interpreted as 
the mean daily geographic longitude of the satellite (or nodal longitude for the near- circular inclined- 
orbit satellite). In a second-order earth-longitude gravity field, Equation 1 is merely a special case 
of Equation 3 for zero inclination. For the second-order case at least, the long term trajectories in 
References 6 and 7 have shown it is proper to replace in Equation 1 by \, the general longitude 
position of the 24-hour satellite in the drift. Thus for second-order gravity drift of an inclined 
24- hour satellite the relevant longitude acceleration regime is 

X = - Tf- |f 22 F(i) 22 sin 2(\-\ 22 )| rad/(sid. day) 2 , (5) 


where 


) 22 


cos 2 i + 1 
2 


( 6 ) 


F( i ) 2 2 is an ’’inclination factor" which may be applied to the zero inclination regime to get the 
proper acceleration regime for the inclined satellite. It is not at all obvious, even for a circular 
satellite orbit, that one can apply such a simple factor, independent of the longitude, to modify the 
equatorial regime. In Reference 6, Equation 6 was found to be a good approximation of the cor- 
rect factor for i < 45° . 

In the first section of the present report it is shown that for the inclined circular orbit of the 
24-hour satellite, the second-order factor, which is exact for all inclinations, is 


F(i) 22 


( 1 + cos i ) 2 
4 


(7) 


It is also shown there, rigorously, that for third and fourth-order as well as second-order longi- 
tude gravity, such exact inclination factors exist to modify the higher order equatorial regimes of 
Equation 1. 

The modification (on an or bit- aver aging basis) gives the complete orbit- aver aged drift re- 
gime of the inclined- orbit satellite (to the fourth order) as 


‘ ^ 17 " [ F nm F ( i )nm sinm (^” X nm)] rad /(sid. day) 2 , (8) 

n = 2 m = 1 
For n— m even 
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where 


g s - "e 2a s 1 ^! a s ■ 

Equation 8 (without the inclination factors) was first proposed as a good approximation for the 
second order drift regime of the equatorial 24- hour satellite by Frick and Garber (Reference 2), 
and for the higher order regimes by R. R. Allan (Reference 4). Earlier than Allan, in Reference 5, 
the author limited his investigation of equatorial 24- hour drift in a fourth- order field to the vicinity 
of the initial longitude placement of the satellite and derived Equation 1 which is seen to be a special 
case of Equation 8. 

DERIVATION OF THE LONG-TERM ORBIT -AVERAGED DRIFT OF A 24-HOUR, NEARLY 
CIRCULAR, INCLINED-ORBIT SATELLITE IN A FOURTH-ORDER EARTH GRAVITY FIELD 

In Reference 6 the technique of orbit averaging the disturbing force along the track of a 24-hour 
satellite to obtain the relevant long-term (secular) drift motion was proposed for the inclined-orbit 
satellite. While the energy -changing disturbing force on the equatorial satellite due to the gravity 
of the earth is constant over an orbit, it is, in general, quite a complicated function of the instantane- 
ous orbit position of the satellite if the orbit plane has the slightest inclination. It can be antici- 
pated, then, that for the inclined- orbit case the exact drift equations of motion for the 24-hour satel- 
lite will contain many small, time-varying "functions" of daily period, introducing non-linearities 
which will be difficult to deal with. Since it is desired at the start to arrive at a theory of long- 
term (i.e., greater than one day) drift motion, it is natural (since the magnitudes of the short period 
terms are small) to seek a perturbation technique which, in a first approximation, smooths out these 
short period effects. The orbit -averaged perturbation force technique chosen in Reference 6 seemed 
a natural one to complement the idea used there of dealing directly with the daily perturbation of the 
two-body energy of the satellite. Computation of the mean daily effect for second-order longitude 
gravity was not too difficult. As a first approximation to the exact drift motion, numerical integra- 
tions over many months proved the long-term theory derived in this manner to be an excellent ap- 
proximation of reality, even in the presence of a more complicated gravity field with much stronger 
effects of daily period. Nevertheless, the technique of direct orbit averaging of perturbation forces 
has its drawbacks. The chief one is computational. It leads immediately to definite integrals (over 
2tt) involving the anomaly 6 of the satellite and the inclination i . (The eccentricity is not involved 
since only the first-order perturbation of a circular orbit is sought.) These are cumbersome to re- 
duce to simple closed forms. Nevertheless, it is easy to show which integrals, no matter how com- 
plicated, will orbit average to zero and which will yield non-zero, secular, long-term effects. The 
cumbersome integrals yielding non-zero mean daily effects can, at worst, be left as they are and 
numerically integrated for any inclination to yield the proper "inclination factor" corresponding to 
that secular effect. It will be found that in every harmonic case (through fourth-order) the inclina- 
tion integral can be directly reduced to yield a factor equivalent to one computed for this harmonic 
by an entirely different method* (Reference 9). Though no necessary equivalence condition has been 
found between the two theories for the development of first-order secular terms of the earth's grav- 
ity disturbing function, it is presumed on the basis of the exact agreement of these through fourth- 
order that such inclination factors will be equivalent in the two theories to all orders of earth gravity. 


*Kaula, W. M., private communication. 
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Proceeding on the basis of the orbit- averaging theory established by Equation 10 in Refer- 
ence 6, the long-term longitudinal drift acceleration of a 24-hour, inclined, nearly circular-orbit 
satellite is given by 


-12 TT 2 F 

7-5- rad/(sid. 

^eA 


day) 2 


( 9 ) 


where 


F = ^ F t ( 0 , i)dfi = F t (&) . ( 10 ) 

f t is the tangential or along-track perturbation force acting on the satellite at an angle G from its 
ascending node. Similarly, the orbit- averaged rate of change of the semi- major axis of the satel- 
lite orbit, coupled with Equation 9, is given by Equation 11 in Reference 6 as 

477 a s 3 F 

a = length units/sid. day . (11) 
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Figure 1 shows the position of the 24- hour 
satellite with respect to the earth and the cel- 
estial sphere. F r , F^ and F x , earth- gravity 
perturbing forces in the radial, latitude, and 
longitude directions, are assumed to be acting 
on the satellite at S (see Appendix A). As long 
as the orbit is nearly circular, F r will make a 
negligible contribution to F T . The tangential 
perturbing force will then be close to the cir- 
cumferential perturbing force and we can write 

F T = F^ cos a + F^ sin a , (12) 


Figure 1— Position of a 24-hour satellite with a near- w h e re F. and F A are the latitude and longitude 

circular orbit with respect to the earth and the celestial , .... 

sphere> components of the gravity force on the satellite 

and a is the angle between a meridian plane and 

the orbit plane when the satellite is 6 distant from its ascending node. Writing out Equation 12 in 

its full fourth- order- gravity form will involve the sum of contributions to F T from every earth 

harmonic (zonal and tesseral) up to and including J 44 (see Appendix A). 

We will now show that the zonal terms of this series (terms in J 20 , J 30 and j 40 ) give zero 
orbit- averaged effects; i.e., F(zonal) = 0. 
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The only force component of zonal gravity capable of having an energy- changing effect on the 
satellite is the latitude one, if the orbit is circular. 

Orbit-Averaged Effect of J 20 

From Appendix A, the relevant term is 


F 0, 20 



{ 3 J20 


sin <£> cos 




(13) 


But from the spherical trigonometry of the orbit on the celestial sphere (Figure 1), the following 
relations must hold: 


sin <$> - sin i sin 9 , 


AL - tan" 1 (cos i cos 6) 
tan 

cosa “ t^9 * 


(14) 


Substitution of Equations 14 and 13 into Equation 12 gives 


F 


T, 20 



tan < fi sin <f> cos eft 
tan 9 


s in 2 4> cos 9 
*20 sin# 



sin 2 i sin 2 6 


(15) 


where k 20 is a constant over the circular orbit assumed. But the integral of Equation 15 over an 
orbit (0<9 < 2 n) is zero, or 


C T,20 


- JL [ 2 

' 2n 

Jo 


Ft, 20 ^ 


= 0 


(16) 


Thus, as long as the orbit remains circular, the J 20 term will produce no change in the energy of 
the satellite over periods of time which are multiples of one day. It must be noted at this point that 
this first-order (zero eccentricity) result for J 20 may not be sufficient in all cases to guarantee 
negligible long-term energy contributions compared to those of longitude gravity for actual orbits 
of small eccentricity. While the effects of small eccentricity can generally be ignored as second- 
order in the theory of the tesseral perturbations, they will be about 1000 times larger in the theory 
to the same order for j 20 due to the relatively large oblateness of the earth compared to the 
other mass anomalies (see DISCUSSION). 
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Orbit-Averaged Effect of J 30 

From Appendix A the relevant term is 


^<£,30 


k 30 sin 2 4> - l) cos , 

where k 30 is a constant for the circular- orbit satellite. As in Equation 15, we now get 


(17) 


F T,30 F <£,30 COS a 


tan 4 > cos <£ , . . 

k 3 0 tan 9 (5 sin 2 <f> - 1) 


k 3Q sin i cos 9 (*> sin 2 i sin 2 0 - l) . 


f 


Since cos 9 d9 - 0, the orbit- averaged tangential force from Equation 18 is 


■ T , 3 0 


5 k 30 sin 3 ; r27T 


2n 


-r 

Jo 


i 2 9 cos 9 d9 - 0 , 


(18) 


(19) 


since sin 2 Q cos 6 is an odd function with respect to 8 - 77/2 and 377/2. 

Orbit-Averaged Effect of i 40 

From Appendix A the relevant term is 


F 0,4O " k 40 sin<£cos [l sin 2 <f> - 3 } 


(20) 


As before, 


F 


T, 40 


4 


But since sin20d£ 


„ ^ \ tan cjb 

F 0 )4 o cosa " k 40 sin<£ cos <£(7 sin J <£ - 3J tan q 

k 40 sin 2 d>cos 6 ^ ^ k 4Q sin 2 i sin 29 

s in Q (7 sin 2 - 3) = 2 (7 sin 2 i sin 2 6 - 3) . 


F 


T, 40 


_1_ 

277 


r2 tt 

F T,40d5 

Jo 


7 k 40 sin 4 i 

477 



sin 2 9 sin 29 d 9 


0 , 


( 21 ) 


because sin 2 6 sin 2 8 is an odd function with respect to 9 = 77 . 
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It is conjectured that similar zero results can be expected to hold for all higher zonal gravity 
forces due to their latitude symmetry. The longitude harmonics of the gravity field of the earth 
will now be treated in a similar manner. 


Orbit-Averaged Effect of j 2 2 

Since geographic longitude dependence is involved in all the longitude effects, it will be neces- 
sary to develop certain trigonometric relationships involving the longitude excursion in the "figure 
8" 24-hour ground track (see Figure 3, Reference 6), as a function of nodal argument 8. 

In Reference 6 (see also Figure 1), it is shown that the excursion A x from the ascending 
equator crossing (at X 0 geographical longitude) in the "figure 8" track is given by: 


AX = AL - co “ tan 1 (cos i tan 6) - 8 


( 22 ) 


Thus from Equation 22 the following identities can be seen to hold: 


cos AX - cos 8 cos [tan 1 (cos i tan + sin 8 sin [tan~ 1 (cos i tan #)] 


cos 2 8 + sin 2 8 cos i 
(cos 2 8 + sin 2 8 cos 2 i ) 1/2 


1 - sin 2 8( 1 - cos i ) 
(1 - sin 2 8 sin 2 i ) 1/2 


(23) 


sin AX - cos 8 sin [tan 1 (cos i tan 8 )J - sin 8 cos [tan " 1 (cos i tan #)J 


cos 8 cos i tan 8 - sin 8 
(l + cos 2 i tan 2 (9 ) 1/2 


sin 20 (co s i ~ 1 ) / 04 ) 

2 (l - sin 2 8 sin 2 i ) 1/2 


cos 2AX - 1 - 2 sin 2 AX - 1 


s in 2 2 8( 1 - cos i ) 2 
2 (l - s in 2 8 s in 2 i) 


(25) 


sin 2 AX - 2 sin AX cos AX 


[l - sin 2 8( 1 - cos i )] [s in 2 #(cos i - 1 )] 
1 - s in 2 8 s in 2 i 


s in 2 <9( 
cos 


cos i - 1) cos i ] sin 2<9(cos i “ 1") 1 /oc\ 

i-rr - 1 * (cosTTTj I t— ^ 7 — ! < 26 > 
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cos 3 AX 


cos(2AX + A X) 


cos 2AX cosAX - sin 2AX sin AX 


■{*- 
1 - sin 2 0(1 - cos i)l 


sin 2 20(1 - cos i) 2 l 

fl - sin 2 6 ( 1 - cos i ) 

2 (l - sin 2 0 sin 2 i)J 

! [ (l ~ sin 2 0 sin 2 i) 1/2 

("sin 20(cos i - 1)1 

j sin 20(cos i - 1 ) 1 

[l - sin 2 6 sin 2 ij 

[2 (l ~ sin 2 0 sin 2 i) 1/2 J 


( 27 ) 


sin 3 AX - sin (2 AX + AX) - sin 2 AX cos AX + cos 2 AX sin AX.. 


[l “ sin 2 0(1 - cos i )] [sin 20(cos i - 1 )] [l - sin 2 0( 1 - cos i )] 
(l - sin 2 6 sin 2 i) (l - sin 2 6 sin 2 i) 1/2 

f sin 2 20(1 "cos i) 2 ] f sin 20(cos i - 1) 1 

\ 2 (l - sin 2 0 sin 2 i) J [2 (l - sin 2 0 sin 2 i) 1/2 J 


cos 4AX - cos 2(2AX) - 1 - 2 sin 2 2AX 


= 1 “ 2 


( sin 20(cos i ~ 1 ) cos i sin 20(cos i - 1 ) V 

t COS 1 + 1 + COS i + 1 ^ _ sin 2 e sin 2 ij 


[l - sin 2 0(1 ~ cos i )] 2 [sin 20(cos i - 1)1 2 

1 - 2 ” r— . . 2 . 12 " ’(28) 

[1 - sin 2 6 sin 2 lj 2 ' * 


sin4AX - sin2(2AX) - 2 sin 2AX cos 2AX 


- 2 {l “ sin 2 0( 1 ~ cos i )} 


sin 20(cos i - 1) 
(l - sin 2 0 sin 2 i) 


^ _ sin 2 20(1 ~ cos i) 2 
2 (l - sin 2 6 sin 


2 0 


(29). 


Continuing with the evaluation of the along- track force due to the longitude and latitude contributions 
from J 22 (see Appendix A), we have 


F T, 2 2 - F A.22 sina + F <M2 COSa 


“ k 22 (sin 2(X. - ^- 22 ) [sin a cos 0] + cos 2 (A. - A. 22 ) [cos a sin 4 > cos 
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where 


k,' = 6 J 2 




122 K ) 2 

When the orbit geometry relations (from Figure 1): 

X — X Q + AX , 


COS 1 

sina - ’ 


( R 0 /O 2 


tan 4 > 

cos a - tST7 ’ 


and 

sin 4> ~ sin i siri 6 

are introduced into this form, the tangential force on the 24-hour satellite due to j 22 is 


(29a) 

(29b) 

(29c) 


(29d) 


F t, 22 k 22 -ssin2(X 0 X 22 ) 


cos i cos 2AX - 


sin 26 sin 2 i sin 2AX 




+ cos 2 (X Q X 2 2 j 


sin 26 sin 2 i cos 2AX 

cos i sin 2AX + - 2 f * (30) 


When the longitude excursion relations, Equations 25 and 26, are used, Equation 30 becomes 


F = k 

r T, 2 2 *22 


_ /. . \ 1 

sin 2 2#(1 - cos i) 2 

sin2(\ 0 -\ 22 W 

* 2 (l “ sin 2 6 sin 2 i) 


sin z 1 


2 [sin 2 20 (cos i - 1 )] 


1 - sin 2 6(1 - cos i) 


1 - sin 2 6 sin 2 i 


r n [sin 2#(cos i - 1)J 

+ cos 2 (;v 0 -\ 22 Wcos i[l- sin 2 0(1 -cos i)] ^ ^ e sin j 7 


sin 2 i sin 26 


sin 2 2 6 ( 1 - cos i) 2 
* 2 (l - sin 2 6 sin 2 i ) 


(31) 
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When the factor of cos 2(\ 0 -\ 22 ) in Equation 31 is examined, it is noted that it is an odd function 
with respect to 6 = tt. Therefore, the orbit- averaged (0 < 6 < 2n) contribution of the cos 2(\ 0 - k 22 ) term 
in the forcing function is zero. It is this result, which will be found to be valid for all the higher 
order tesserals, which enables one to write the relevant long-term drift equations for the inclined- 
orbit 24-hour satellite as a simple (inclination- factored) modification of the equations for the 
equatorial satellite. 

Thus, only the factor of sin2(k 0 ~k 22 ) in Equation 31 need be considered, or 

F Ti22 3 k 22 sin 2(\ 0 - k 22 ) F(0, i) 22 , ( 31a ) 


where 


F(0, i) 22 _ cos i ~ sin 2 2 6 


cos i - sin 2 2 6 


( 1 ~ cos i ) 2 cos i + sin 2 i(cos i - 1 ) + sin 2 i sin 2 6(1 - cos i ) 2 
2 (l - sin 2 6 sin 2 i) 

( 1 - cos i ) [cos i( 1 - cos i ) - sin 2 i + sin 2 i sin 2 6(1 - cos i )] 
2 (l - sin 2 6 sin 2 i ) 


cos i - sin 2 20(1 - cos i) 


cos i ~ cos 2 i - sin 2 i + sin 2 i sin 2 0(1- cos i ) 
2^1 - sin 2 6 sin 2 i) 


* + sin 2 2 6( 1 - c os i) 2 [l si n i s in o\ _ cos i + sin 2 2 6( 1 ~ cos i) 2 
2 1 - sin 2 i sin 2 6 2 


(32) 


Thus, the orbit-averaged along- track perturbation force due to J 22 is (from Equations 31a and 
32) 


Fj, 2 2 ^22 s i n 2(Aw Q ^- 22 ) F(0» i) 22 


[ k 2 2 Sin2 (^0 ^ 22 )] 

1 f 2n f 

2^7 I cos i + (1 - cos i) 2 

k 22 sin 2 (k Q - \ 2 2 ) 

f . (1 - cos i ) 2 1 

[cos i + ^ 4 J 

sin 2 (k Q - k 2 2 ) 
k 22 4 

[4 cos i + 1 - 2 cos i + cos 2 i] 

sin 2 (k 0 - k 2 2 ) 
k 22 4 

[l+2 cos i + cos 2 i] 

sin 2 (\ 0 - ^ 22 ) 

- k 22 4 (1 + cos 1) 2 


20 ] 


d 6 


( 33 ) 
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It is seen that the orbit- averaged, energy- changing perturbation force on an inclined-orbit 24-hour 
satellite due to J 22 is just F(i) 22 times that on the equatorial satellite (Equations 1 and 2, or A-4) 
at the same mean longitude K 0 , where (from Equation 33) 


F Ci ) 22 


( 1 + cos i ) 2 


(34) 


Orbit-Averaged Effect of J 3) 

The tangential force on the 24-hour inclined- orbit satellite at Q from its ascending node, due 
to J 31 , is (see Appendix A) 


1 T,3 1 


( 6 ) = 


F A,3l Sina + 


0,3 1 


n ^sin (A. - \ 31 ) [5sin 2 c£-l] sin a + cos (A - A 3 [l5 sin 2 - 1 1] sinc^cosaj- 


(35) 


where 


k 


3 1 


I ( R o/a s ) 3 J31 • 


(35a) 


Introduction of the orbit geometry relations, Equations 29a through 29d, causes Equation 35 to 
become 


F T, 3 1 R 3 1 "A S ^ n (^0 ^3l) 


(5 sin 2 6 sin 2 i - l) cos i cos AA 


(l - sin 2 i sin 2 6) 1/2 


(15 sin 2 i sin 2 9 - ll) (sin 2 i sin 2d) sin AA 


(l - sin 2 i sin 2 &) 1/2 


+ cos(X. 0 -\ 31 ) 


(5 sin 2 6 sin 2 i - l) cos i sinAA 

(I- 


sin 2 i sin 2 6?) 1/2 


(l5 sin 2 i sin 2 Q ~ ll) (sin 2 i sin 2 &) cosAA 
2 (l " sin 2 i sin 2 9j 1/2 


( 36 ) 
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Introducing the expressions for the longitude excursion in the "figure 8", Equations 23 and 
24, causes Equation 36 to become 


F T,31 (0) k 31 "S sin (^0 ^3l) 


cos i {l — sin 2 8 [1 - cos i]} (5 sin 2 8 sin 2 i - l) 


1 - sin 2 i sin 2 8 


sin 2 i sin 2 2#(cos i - 1) (l5 sin 2 8 sin 2 i ~ ll) 

4 (l - sin 2 8 sin 2 i) 


+ COS (\„-^3l) 


cos i (5 sin 2 8 sin 2 i ~ l) sin 2#(cos i - 1 ) 
2 (l - sin 2 8 sin 2 i) 


sin 2 i sin 20 (l5 sin 2 i sin 2 8 - ll) {l - sin 2 8 { 1 - cos i )} "1 
2 fl - sin 2 8 sin 2 i) j 


Once again it is seen that the factor of cos (\ 0 - a . 31 ) in Equation 37 is an odd function about 8 = n 
so that its orbit- averaged effect over 0 < Q < 2n is zero. When only the relevant term sin (\ 0 - \ 31 ) 
term is retained, Equation 37 becomes 

k 3 1 s xn A.3 j j ^ 

F T,3l( 0 ) = 4(l- sin 2 8 sin 2 i) \ 4 cos 1 t 1 " si ° 2 *(l-cosi)] [5 sin 2 8 sin 2 i - l] 


sin 2 i sin 2 20(cos i - 1) [l5 sin 2 8 sin 2 i -ll]^ (37a) 


k 3 j s in ^-3 J cos i r -j r 

4(l - sin 2 6 sin 2 i ) 1 1 + cos i [ cos 1 + (* " si ° 2 6 si ° 2 OJ [“ ( X " sin2 6 sin2 £ ) + 4 si " 2 6 sin2 2 J 


sin 2 i sin 2 26(cos i - 1) [-10 (l - sin 2 i sin 2 d) + 5 sin 2 


1] \ (37b) 


J-4 cos i (l - sin 2 8 sin 2 i) -4 cos 2i 16 cos i sin 2 8 sin 2 i 10 sin 2 i sin 2 2#(cos i 
k 3 i sin (\ Q _ ^3i)S 4(1 + cos i) 4(1 + cos i) + 4(1 + cos i) + 4 


16 cos 2 i sin 2 i sin 2 8 ' * _ n 1 

l + cos i - sin 2 i sin 2 20(cos i ~ 1) [5 sin 2 8 sin 2 i — 1 J I (37c) 

4 (l - sin 2 8 sin 2 i) | 


[ 16 cos 2 i sin 2 i sin 2 8 + sin 4 i sin 2 28 (5 sin 2 8 sin 2 i - l) 

k 3 2 sin (/V Q \ 3 j) <A(6) + 4(1+ cos i ) (l - sin 2 8 sin 2 i) 
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- k 31 sin \ 0 -\ 31 <A(0) 


16 cos 2 i sin 2 i sin 2 (9 - sin 4 i sin 2 26 j^(l - s in 2 6 sin 2 i) - 4 sin 2 6 sin 2 i| 


4(1+ cos i ) (l - sin 2 6 sin 2 i) 


>(37e) 


B(6) 


k 31 sin ( K 0~ K 3l) 


, r . '■ 

I “4 cos i (1 - sin 2 6 sin 2 i) - 4 cos 2 i + 16 cos i sin 2 6 sin 2 i - 11 sin 4 i sin 2 26 

4(1+ cos i ) 

_ J 

V 


16 cos 2 i sin 2 i sin 2 (9 + 4 sin 6 i sin 2 2 6 sin 2 6 

+ -7 7 y (37f) 


4(1+ cos i ) (l - sin 2 6 sin 2 i) 


sin 2 6 sin 2 i U (l - sin 2 i) + sin 4 i sin 2 26 


(37g) 


k 3 1 sin (\ 0 -k 3l ) <B (6) + 


sin 2 6 sin 2 i ^4 - 4 sin 2 i (l - sin 2 i sin 2 6) - 4 sin 4 6 sin 4 i 


( 1 + cos i ) (l - sin 2 0 sin 2 i) 


(37h) 


C (&) 


4 sin 2 9 sin 4 i 4 sin 2 0 sin 2 i [l - sin 4 8 sin 4 i] 

k 31 Sln (\> “N>l) 1“ COS i 






( 1 + cos i ) (l - sin 2 6 sin 2 i) 




4 sin 2 6 sin 2 i (l + sin 2 6 sin 2 i) 
k 31 Sin ( X 0 ~^3l) i C (^) + 1 + COS i 


-J 


(37i) 


(37j) 


Reducing D(<9) in Equation 37j gives 


D(^) 


4 sin 2 6 sin 2 i(4 cosi + cos i) - 4 cos i(l + cos i) - 16 sin 2 6 sin 4 i - 11 sin 4 i sin 2 26 + 16sin 2 £sin 2 i (l + sin 2 6 sin 2 i) 

4(1 +cos i ) 


sin 2 6 s in 2 i(5 cos i + 4) 
( 1 + cos i ) 


sin 4 i 

cos i - 4 ( 1 + cos i) {ll sin 2 2 6 - 16 sin 4 6 + 16 sin 2 


(37k) 
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From Equations 37j and 37k the orbit-averaged tangential perturbing force on the inclined- orbit 
24- hour satellite, due to J 31 is 



f sin 2 i(5 cos i + 4) 

= k 3 i sin (X 0 -^ 3 i) V cos 1 + 2(1 + cos i) 


sin 4 i 
4(1+ cos i ) 



- k 31 sin 


(^0 ^ 3 l ) 1 


(1 - cos i) (5 cos i +4) 15(1 ~ cos i) sin 2 


/ \ 5 cos i 5 /„ 

- k 3] sin (\ 0 -k 31 ) <- cos i + 2 ' 2 l 1 sln 1 

1 F(i ) 31 

A 


« 15 . , . 15 . 2 . 

) + 2 " 2 cos l - "g" sin" 5 i + g sin 1 cos 1 


N - 

x |l+cosi 5 sin 2 i(l + 3 cos i) 

-k 3 ,sin (A. 0 -*-3lM 2 8 

U J- 


(38) 




But -k 31 Sin (\ 0 -\ 31 ) = F 31 sin2(>v 0 -\ 31 ) , which is the tangential perturbation force on the 24- 
hour equatorial satellite at longitude k 0 . Thus, the orbit-averaged force on a 24-hour inclined 
circular orbit satellite whose ascending node is at is just F(i) 31 times the force on the 
equatorial 24-hour satellite at ^ 0 , due to J 31 , where 


F( i >3 1 


1 + cos i 
2 


5 sin 2 i 
8 


(1 + 3 cos i ) 


(39) 


Orbit-Averaged Effect of J 32 

The tangential force on the 24-hour inclined- orbit satellite at Q from its ascending node, due 
to J 32 , is (see Appendix A) 


F T ,32 ( ^ “ F A,32 Sina + F ^,32 COSa 


k J2 -|sin 2 (>t - ^ 32 ) [2 sin a sin<£ cos + cos 2 (A. - \ 32 ) |cos a (3 sin 2 l) cos j. 


(40) 


where 


3 2 


= 15 


( 3 s) : 


( R o /-,) 3 J 


32 
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Introducing the orbit geometry relations, Equations 29a through 29d, causes Equation 40 to become 


F t 32 (0) - k 32 |sin 2 (^ 0 - A 32 ) |^2 cos i sin i sin 0 cos 2 A K - cos 0 sin i (3 sin 2 i sin 2 9 - l) sin 2 AxJ 

+ cos 2 (^o ^ 32 ) [ 2 cos i sin i sin 0 sin 2A\ + cos 1 9 sin i (3 sin 2 i sin 2 9 - l) cos 2A\jJ > *(41) 

When Equations 25 and 26 are used, Equation 41 now becomes 


1 T, 3 2 


r 



sin 2 (\„ -\ 32 ) 

sin 2i sin 0 « 

sin 2 20( 1 “ cos i ) 2> | 

^ 2 (l - sin 2 0 sin 2 i)J 

l 




6 (3 sin 2 i sin 2 9 - l) 


{l - sin 2 0(1 - cos i)} {sin 20 (cos i - 1 )} 


1 - s in 2 9 sin 2 i 


+ cos 2 (K Q A -3 2 ) 


sin 2i sin 9 


{l - s in 2 0 ( 1 - cos i )} {s in 20 (cos i - 1 )} 


1 - sin 2 9 sin 2 i 


] sin 2 20(1 ~ cos i ) 2 

+ cos 0 sin i (3 sin 2 i sin 2 0 - l) "S 1 * ~T ‘ 

v ' 2 (l — sin 2 0 sin 2 1 ) 


(42) 


The factor of sin 2(a 0 - \ 32 ) in Equation 42 is an odd function with respect to 0 = n. Thus it has a 
zero orbit- averaged effect. The factor of cos 2 (\ 0 - \ 32 ) in Equation 42 is an odd function with 
respect to 0 = 77/2 and 3tt/2. The orbit- averaged effect of this term is also zero. Therefore, the 
orbit- averaged contribution to the longitude acceleration of the inclined- orbit 24-hour satellite due 
to the tesseral J 32 is zero, as it is at all orbit anomalies, identically, for the equatorial satellite. 


Orbit-Averaged Effect of i 33 

The tangential force on the 24-hour inclined- orbit satellite at 0 from its ascending node, due 
to J 33 is (see Appendix A) 


- F x ,, sina + F , , , cos a 


r T,33 r A. , 3 3 


0,3 3 


= k 33 jsin3(\-\ 33 ) sin a cos 2 <t> + cos 3 (A. - X> 33 ) cos acos 2 <£sin 


(43) 


where 




45 (a ) 2 ( R °/ as ) 3 ^33 ’ 
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Introducing the orbit geometry relations, Equations 29a through 29d, causes Equation 43 to become 
F T 33 (#) = k 33 -|sin 3 (\ 0 - A.33) |cos i (l - sin 2 8 sin 2 i ) 1/2 cos 3AA. - — “ sin 20 (l - sin 2 0 sin 2 i) 1/2 sin 3A^j 


+ COS 


3 (X 0 “ ^33) [cos i (l - sin 2 0 sin 2 i) 1/2 sin 3A\ + “ " sin 20 (l - sin 2 6 sin 2 i) 1/2 cos 3A\|^(44) 


Using the longitude excursion relations, Equations 27 and 27a, in Equation 44 gives 


^T,33 ^33 3(^o ^33) 


COS i <1 


sin 2 20 ( 1 - cos i ) 2 


(T^esin*i)J {l “ sin 2 0(1 ~ cos i)} 


f sin 2 2#(cos i - l) 2 

cos i (l - s in 2 0 ( 1 - cos i )\ < — 7 „ ~ — r 

1 ;; 1 2 (l -sin 2 £sin 2 i) 


\26 


[l - sin 2 0(1 ~ cos i)j [sin 2#(cos i - 1)] [l ~ sin 2 0(1 ~ cos i)] 
(l - sin 2 0 sin 2 i) 


-j2 £ s i n 20 f sin2 20(1 " cos i) 2> [ (sin 20(cos i - 1) 
' * 2 (l - sin 2 0 sin 2 i) 


. + cos 3(\ 0 -\ 33 ) 


J"[l - sin 2 0(1 - cos i ) j [sin 2 < 9 (cos i — 1 )] [l - sin 2 0(1 - cos i )] 
(l - sin 2 0 sin 2 i) 


+ cos i il ~ 


sin 2 20 ( 1 ~ cos i ) 2 [ [sin 2 ( 9 (cos i - 1 ) 


2 [l “ s in 2 0 s in 2 i) 


sin 2 i sin 2(9 f sin2 20(1- cos i) 2 

+ o -<1 


° si)2 l r , 

2 (l - sin 2 0 sin 2 i) J U ~ ^ 0(1 - cos i)} 

sin 2 i sin 20 J [* “ sin2 d ( l ~ cos *•)] [ sin 20(cos i - 1)] [sin 26 ’(cos i - 1)] 

2 (l - sin 2 0 sin 2 i) 


' (45) 


Once again the factor of cos 3 (^ 0 -^ 33 ) above is an odd function with respect to 0 = v, so that it 
orbit-averages to zero. 
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Rewriting Equation 44 and including only the relevant sin 3 (k 0 - k 33 ) term gives 


F T ,33 k 33 Sin3 (^0 ^ 33 ) ^ cos 1 


sin 2 20 ( 1 - cos i ) 2 

2 (l ~ sin 2 6 sin 2 i )J 


[l - sin 2 6(1 - cos i )] 


[sin 2 2#(cos i - l) 2 } s i n 2 i [l “ sin 2 6(1 ~ cos i)] 2 sin 2 2 < 9 (cos i 

cos ill- sin 2 6(1 - cos i )1 — 7 — — r~ - — ^ , "* 

1 -* 2 (l - sin 2 6 sin 2 i) 2 (l - sin 2 0 sin 2 i) 


4 1 [sin 2 2#(cos i - 1)] 


1 - 


sin 2 20(1 ~ cos i) 2 
2 (l “ sin 2 6 sin 2 i ) J 


\ (*) 


k 33 sin 3 (\ 0 ~ ^33) < cos i [l - sin 2 0 ( 1 - cos i )] - 4 [sin 2 20(cos i - 1 )] 

— — J 


sin 2 20 ( 1 - cos i ) 2 


cos i {l - sin 2 6(1 - cos i )} + 2 -[sin 2 2 < 9 (cos i - l)j 


2 (l - sin 2 6 sin 2 i) 

- cos i {l - sin 2 6(1 ~ cos i)} ” (cos i - 1) {l “ sin 2 6(1 - cos i)} 
sin 2 20 ( 1 ~ cos i ) 2 


^ 33 S ^ n ^(^0 ^"3 3 ) 1 ^1 ^ 


2 (l - sin 2 6 sin 2 i) 


-2 cos i 

(1 +cos i) { cos 1 + 1 “ sin2 6 sin 2 i] 


sin 2 i sin 2 2#(cos i - 1) [cos i + (l - sin 2 6 sin 2 i)} 2 
4 + 1 + cos i 


f sin 2 20(1 - cos i) 2 

k,, sin 3 (A. n -A.,.) < A. (0) + 7 — < 

33 v 0 33 ' | i v ’ 2(l- sin 2 0 sin 2 ij 


cos i + (l - sin 2 0 sin 2 i) 

1 + c^TI f 2 cos i 




— 1 

'i 

|cos i + 

L ~ sin 2 6 sin 2 i)J 

) sin 2 i „ 

> + 4 sin 2 2<9(cos i - 1 ) 

l 

J 


r 


B^(0) 


sin 2 20( 1 - cos i ) 2 


k 33 sin3(\ 0 *- 33 ) ^ A i ( 0 ) + 2(1+ cos i ) [l - cos i - sin 2 0 sin 2 i] + 


sin 2 20(1 




2(1 


1 + cos i cos 2 + (l “ sin 2 ® sin 2 i)} + sin 2 i sin 2 0 (l - sin 2 0) (1 - cos i) 


1) 


(46) 


(46a) 


(46b) 


(46c) 


- cos i 
0 sin 2 

(46d) 
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{{ 


c, ( e ) 


} sin 2 20(1 - cos i) 2 cos i sin 2 20(1 - cos i) 2 

k 33 sin3 (x. 0 — X 33 ) (g)+ 2(1 + cos i ) J + 2(1 +cos i) (1 -sin 2 0sin 2 i) 




cos 2 i - sin 4 i sin 2 8 (l - sin 2 #)J^- (46e) 


sin 2 26(1 - cos i) 2 
k 33 sin 3(^ 0 ^33) 'jCj (0) 2(1 + cos i ) (l - sin 2 6 sin 2 i) 


|~1 - (l - sin 2 6 sin 2 i) sin 2 i - sin 4 i sin 4 j* (46f) 


Dj (0) 


sin 2 2^(1 - cos i) 2 (l + sin 2 6 sin 2 i) sin 2 i sin 2 20(1 - cos i) 2 l . 
k 33 sin3(X- 0 “^ 33 ) i c x (&) ~ “ 2( 1 + cos i ) + 2(l+cosi) ' g 


Thus, the orbit- averaged tangential force due to j 33 is 


. r 1 C 27T 

33 (^) = ~2rr F T,33^) d ^ = k 33 sin 3 (^o “ ^ 33 ) 2tt ] ( 0 ) d 0 

’ Jo 


/ X J r ( 1 “ COS i )~| sin 2 i 

= k 33 sin 3(^ 0 -X 33 ) |cos 1 1^1 2 J 8 


(1 - cos i ) 3 
(cos i“l) + 4(1+ cos i ) 


(1 - cos i) 2 sin 2 i (1 - cos i) 2 cos i ( 1 - cos i) 2 
8(1 + cos i ) + 4(1 + cos i) ~4(l+cosi) 


sin 2 i( 1 - cos i ) 2 sin 2 i( 1 - cos i )“ 


= k 33 sin 3(\ 0 -^ 33 ) -\~2~ ( 1 + cos 1) + 


8(1 + cos i ) 4 

4(1 + cos i ) 

(1 - cos i) 2 (1 + cos i) 


) + 8 


1 -j (1 - cos i ) 2 cos i 

( 1 - cos i ) 2 

[4( 1 “ cos i )J + 4(1+ cos i ) 

4(1 + cos i ) [ 


(47) 


(47a) 
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(47b) 


, % j isv si (1 ~ cos i) 2 

= k 33 sin3(\ 0 -\ 33 ) y ~ 2 (1 + cos i) + g(l + cos i) [( 1 + cos i ) 2 + 2(l~cos i) - 2 + 2 cos i] 


v j cos i( 1 + cos i ) ( 1 - cos i ) 2 ( 1 + cos i ) 2 

k 33 sin 3(^0 “^33) 1 2 + 8(1+ cos i) 


k 3 3 sin 3 


( K o ~ ^ 3 3 ) [4 cos i + (1 - cos i) 2 ]j. 


k 3 3 s i n 3 (^0 - ^ 3 3 ) 1 8 [4 cos i + 1 - 2 cos i + cos 2 i] 


, 4 | ( 1 + cos i ) 3 

k 33 sin3 ( X 0' X 33) 1 8 


(47c) 


But k 33 sin 3(^ 0 -\ 33 ) = F 31 sin3(A. 0 -\ 33 ) , which is the tangential perturbation force on the 24-hour 
equatorial satellite at longitude k 0 . Thus, the orbit- averaged perturbing force on a 24-hour, 
inclined- circular orbit satellite whose ascending node is at k 0 is just F(i) 33 times the perturbing 
force on the equatorial 24- hour satellite at k Q9 due to J 33 , where 


F<i > 3 


( 1 + cos i ) 3 
“ 8 


(48) 


Orbit-Averaged Effect of J 41 

The tangential force on the 24-hour inclined- orbit satellite at Q from its ascending node, due 
to J 41 , is (see Appendix A) 


(*) 


r A,4l 


sin a + F 


<t>,4 1 


k 4 1 ‘jsiH “ ^ 44 ) sin a [7 sin 2 cf> ~ 3] sin<£ + cos (A. - ^ 41 ) cos a [28 sin 4 <f> - 27 sin 2 + 3]}, (49) 


where 


k 


41 


s ^ 

* W 
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Introducing the orbit geometry relations, Equations 29a through 29d, causes Equation 49 to take 
the form 


F T,4I (0) " k 40 sin K“^4l) 


cos i (7 sin 2 0 sin 2 i - 3 ) sin 0 sin i cos A A. 


(l - sin 2 0 sin 2 i) 1/2 


sin i cos 6 (28 sin 4 0 sin 4 i - 27 sin 2 0 sin 2 i + 3) sin A k 
(l - sin 2 6 sin 2 )i 1/2 


cos (k Q - \ 41 ) 


cos i (7 sin 2 0 s in 2 i - 3 ) sin 0 sin i sin AA. 
(l - sin 2 6 sin 2 i) 1/2 


(28 sin 4 0 sin 4 i - 27 sin 2 6 sin 2 i + 3) 

+ sin i cos 6 ~ y _ „ . cos AA. 

(l - sin 2 6 s in 2 i) 1/2 


(50) 


Next the longitude excursion relations, Equations 23 and 24, in Equation 50 give 


F T,4 1 k 4 1 S ^ n (^0 ^4l) 


cos i (7 sin 2 6 sin 2 i - 3 ) sin 6 sin i {l - sin 2 0(1- cos i )} 
(l - sin 2 6 sin 2 i) 


sin i cos 6 (28 sin 4 0 sin 4 i - 27 sin 2 6 sin 2 i + 3) sin 20(cos i - 1) 
2 (l - sin 2 6 sin 2 i) 


+ COS (\ 0 -\ 41 ) 


cos i (7 sin 2 6 sin 2 i - 3) sin 8 sin i sin 20(cos i - 1) 
2 (l “ sin 2 6 s in 2 i) 


sin i cos 0(28 sin 4 9 sin 4 i - 27 sin 2 9 sin 2 i + 3) 

+ - . 

(l - sin 2 8 sin 2 i) 

The first term of the factor of sin (a. 0 - A. 41 ) in Equation 51 is an odd function about 8 - n. Thus, 
this term orbit averages to zero. The second term of the factor of sin (\ 0 - X. 41 ) in Equation 51 
also is an odd function about 8 - n as it is controlled by sin 28 cos 9 = 2 sin 8 cos 2 0, which is odd 
about 0 = 77. The first term of the factor of cos (\ 0 -\ 41 ) in Equation 51 is an odd function about 
0 = tt/2 and 377/2, as it is controlled by sin 0 sin 20 = 2 sin 2 0 cos 0, which is odd about 0 = 77/2 and 
377/2. Therefore, this term orbit averages to zero. The second term of the factor of cos (\ 0 -\ 41 ) 
in Equation 51 is also an odd function about 0 = 77/2 and 3n/2. Thus, the orbit average of the entire 
forcing function F T 41 (0) is zero; 

l f 27T 

F T, 4 1 (^) 277 I F T 41 (^) d0 - 0 . 

Jo 


{l - sin 2 0(1- cos i )} 


(51) 
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Thus, there can be no long-term contribution to the longitude drift of a 24-hour, inclined, circular- 
orbit satellite due to this gravity harmonic. 

Orbit-Averaged Effect of J 42 

The tangential force on the 24-hour inclined- orbit satellite at 6 from its ascending node, due 
to J 42 , is (see Appendix A) 


F r ( 42 < e ) ~ F A.42 sina + F <£,42 COS a 


k 42 -jsin 2(/V - \ 42 ) sin a (7 sin 2 4> ~ l) cos <f> 

+ cos 2 (a - k 4 2 ) cos a (7 s in 2 <f> ~ 4) 2 cos c^sin 


(52) 


where 




^42 ( a ) 2 ^ K °^ aS ^ ^ 42 * 

Introducing the orbit geometry relations, Equations 29a through 29d, causes Equation 52 to become 
f t 42 (0) - k 42 jsin 2 (^ 0 - ^ 42 ) [cos i (7 sin 2 6 sin 2 i - l) cos 2AA - sin 2 i sin 2 26 (7 sin 2 6 sin 2 i - 4) sin 2A\J 

+ cos 2 (k Q - \ 42 ) jjcos i (7 sin 2 6 s in 2 i - l) sin 2AA. 

i sin 26 (7 sin 2 6 sin 2 i - 4) cos 2 A\J^ • (53) 


+ s in^ 


Using the longitude excursion relations, Equations 25 and 26, makes Equation 52 become 


F T, 4 2 ^42 V in2 ( X 0 ^ 42 ) 


sin 2 26(1 ~ cos i) 2 

cos i (7 sin 2 6 sin 2 i ~ l) Si ” ~T . 2 Q . 2 .\ 
' ' 1 2 (1 - sin 2 9 sin 2 1 ) 


- sin 2 i s in 26 {l sin 2 6 sin 2 i - 4) 
+ cos 2(k 0 -k 42 ) 


-[l - sin 2 6(1 ~ cos i)} {sin 2<9(cos i - 1)} 
(l “ sin 2 6 sin 2 i) 


cos i (7 sin 2 6 sin 2 i - l) 

T . 2 . 2 {l - sin 2 6(1 - cos i)}{sin 20(cos i - 1)} 

(1 - s in 2 6 sin 2 ij 1 ; 1 J 


sin 2 26(1 ” cos i) 2 

sin 2 i sin 26 (' 7 sin 2 6 sin 2 i - 4) < 1 - f ~ 7~7\ 

\ 2(1 ~ sin 2 6 sin 2 l) 


(54) 
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The factor of cos 2(\ 0 - X 42 ) in Equation 54 is an odd function about 9 = it. Thus, it orbit averages 
to zero and has no relevance on the long term drift of the satellite. 

Considering only the relevant sin 2(\ 0 -^ 42 ) term, Equation 54 becomes 


a 2 ( 0 ) 



r 


A 

Ft.42 <«> 

= k 42 <sin2(A. 0 -\ 42 ) 

cos 

i (7 sin 2 9 sin 2 i - l) 


l 


V 

cos i (7 

sin 2 9 sin 2 i - l) sin 2 2 9 (l - cos i) 

2 - 2 sin 2 i sin 2 29 (cos i 


2 (l _ s in 2 9 sin 2 i) 


(55) 


- k 42 sin 2(^ 0 -A. 42 ) A 2 ( 0 ) 


+ k 42 sin2(\ 0 -\ 42 ) 


cos i (7 sin 2 0 sin 2 i _ l) sin 2 20(1 - cos i ) 2 
2 (l “ s in 2 9 sin 2 i) 


2 sin 2 i sin 2 2 9( 1 - cos i ) 2 ^7 sin 2 9 sin 2 i - l) - 3 ] [l - sin 2 9(1 - cos i )] 
2(l - sin 2 £ s in 2 i) (1 - cos i) 


(55a) 


- k 42 sin 2 (A. q -\ 2 ) A 2 (9) 

k 42 sin 2(^ 0 - \ [ f 9 sin 2 i 

<sin 2 29 ( 1 - cos i ) 2 (7 sin 2 9 sin 2 i - l) - cos i + M - cos j \ {l “ sin 2 9( 1 - cos i )} 

2 (l “ sin 2 9 sin 2 i)l L 


6 sin 2 i s in 2 2 9( 1 - cos i ) [l - sin 2 9( 1 - cos i )] > (55b) 


B 2 ( 0 ) 




- k 42 sin 2 (A- 0 - A 42 ) A 2 (9) + k 42 sin 2 (A . 0 “ A 4 2 ) <sin 2 2 9( 1 - cos i) 2 (7 sin 2 9 sin 2 i - l) 




cos i sin 2 29 ( 1 - cos i ) 2 (7 sin 2 9 sin 2 i — l) - 6 s in 2 i sin 2 29 ( 1 - cos i ) [l - sin 2 9(1 - cos i)] 

2 (l ~ sin 2 9 sin 2 i) 


(55c) 


- k 42 sin 2 (A q -\ 42 ) [A 2 (0) + B 2 (0)] 

k 42 sin 2 (K q - \ 42 ) sin 2 29(1 - cos i) 2 


2 (l - sin 2 9 sin 2 i) 


6 sin 2 i [cos i + (l - sin 2 9 sin 2 i)J I 
cos i (7 sin 2 0 sin* i - l)- ;wl + f(55d) 


( 1 - cos i ) ( 1 + cos i ) 
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= k 42 sin 2 (\ 0 -\ 42 )[A 2 (0)+B 2 (0)] 


k 42 sin 2(\ 0 - ^ 42 ) sin 2 20(1 - cos i ) 2 | . .* , ...... .1 

2 (1 - sin 2 0 sin 2 i) 1 L" ^ " S1 " 9 S1 " + 6 510 9 S1 " X J 


- 6 jcos i + (l - sin 2 0 sin 2 i)jj> (55e) 


= k 42 sin 2 (\ 0 -X 42 )[A 2 (0) +B 2 (0)] 

k 42 sin 2 (A-q - ^ 42 ) sin 2 2&( 1 - cos i ) 2 
2 (l - sin 2 6 sin 2 i) 


<^[l - sin 2 6 sin 2 i] [(- 6 “ cos i) - 6 cos i]j- 


C 2 (0) 




f 'siin 2 2 1 - cos i ) 2 [6 + 7 cos i 

k 42 sin 2(^o -\ 2 )^ A 2 (^) + b 2 (^) " - 2 1 

L ^ - 


"V 


Thus, the orbit- averaged tangential force due to J 42 is (from Equation 55g) 


(55f) 


(55g) 


F t, 4 2 ( 0 ) = ^J o F T , 42 (0)d0 = k 42 sin2(\ 0 -\ 42 ) (^)J o [A 2 (0) + B 2 (0) + C 2 (0>] 


dd 



= k 42 sin 2 (^ Q - \ 42 ) 2t7cos i +7Tcos i (7 sin 2 i) - tt( 1 - cos i ) 2 


7 sin 2 i( 1 - cos i ) 2 J s in 2 20 sin 2 0 d0 - (1 - cos i) 2 (6 + 7 cos i)J- 


(56) 


k 42 sin 2 




(1 - cos i ) 2 


7 sin 2 i , ( 1 -cosi ) 2 

+ 4 ( 1 - cos i) - 4 (6 + 7 cos i ) 


(56a) 


k 42 sin 2(\ 0 \ 42 ) 


■|l4 cos i (l - cos 2 i) - 4 cos i - 2 (l ~ 2 cos i + cos 2 i) 


+ 7 (l _ cos 2 i) (l ~ 2 cos i + cos 2 i) - (l - 2 cos i + cos 2 i) (6 + 7 cos i )j- (56b) 
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k 42 sin 2 (^ 0 -\ 42 ) 

4 {“ 14 cos i + 14 cos 3 i + 4 cos i + 2 - 4 cos i + 2 cos 2 i - 7 


+ 14 cos i - 7 cos 2 i + 7 cos 2 i “ 14 cos 3 i + 7 cos 4 i + 6 ~ 12 cos i 


+ 6 cos 2 i + 7 cos i - 14 cos 2 i + 7 cos 3 i}(56c) 


k 42 sin 2 (X 0 “ ^ 42 ) 

4 {l - 5 cos i- 6 cos 2 i + 7 cos 3 i + cos 4 i} 


(56d) 


k 42 sin 2(^ 0 -k 43 ) 


-jl + 2 cos i + cos 2 i - 7 cos i (l + cos i - cos 2 i - cos 3 i)) 


(56e) 


' j ( 1 + cos i ) 2 7 cos i (l - cos 2 i) ( 1 + cos i )^| 

k 42 sin2(\ -\ 2 )S 4 ~~ 4 I 


(56f) 


1(1 + cos i) 2 7 sin 2 i cos i( 1 + cos i ) 

k 4 2 sin 2 (X 0 -\ 2 )< 4 - 4 


(56g) 


But - k 42 sin 2(\ 0 - \ 42 ) = F 42 sin 2 (\ 0 - \ 42 ) is the tangential perturbation force on the 24-hour 
equatorial satellite at longitude k Q . Thus, the orbit-averaged perturbing force on a 24-hour, 
inclined circular orbit satellite whose ascending node is at \ 0 is just F(i ) 42 times the perturbing 
force on the equatorial 24-hour satellite at k 0 , due to J 42 , where 


F(i) 4 


(1+cosi) 2 7 sin 2 i cos i( 1 + cos i ) 


(57) 


Orbit-Averaged Effect of J43 

The tangential force on the 24- hour inclined- orbit satellite at 6 from its ascending node, due 
to j 43 , is (see Appendix A) 


F Tj 43 (0) ‘ F A,43 Sina + F <M3 C ° Sa 


f 


= k 43 < sin 3 (\ - A 43 ) sin a cos 2 0sin <56 + cos 3 (k - A 43 ) 


cos a (4 sin 2 </> - l) cos 2 4> 


(58) 
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where 


k 43 = 315 ^J 2 (R 0 a s ) 4 J 43 . 

Introducing the orbit geometry relations, Equations 29a through 29d, causes Equation 58 to become 


F T,43 ( 0 ) 


k 43 jsin 3 (^ 0 - ^ 43 ) jcos i (l - sin 2 0 sin 2 i ) 1/2 sin 0 sin i cos 3A K 


S " 3 ° S (4 sin 2 0 sin 2 i - l) (l - sin 2 0 sin 2 i ) 1/2 sin 3 A\J 
cos 3 (k Q - ^ 43 ) Jcos i (l " sin 2 0 sin 2 i ) 1/2 sin 0 sin i sin 3AX. 


+ S_ (4 sin 2 0 sin 2 i - l) (l - sin 2 0 sin 2 i ) 1/2 cos 3A\ 


Using the longitude excursion relations, Equations 27 and 27a, causes Equation 58 to become 


F T,43 ( 0 ) ^43 ^43) 


cos 1 sin 1 s 


in i sin 0 


sin 2 20(1 - cos i ) 2 I 

1 “ -on ' . 2 *\ r {l “ s in 2 0 ( 1 - cos i )} 

2 (1 - sin 2 0 sin 2 1 ) 1 


f ^ |sin 2 20 (cos i - l ) 2 

i sin i sin 0 {l - sin 2 0(1 - cos i)}< 7 " . ' . . \ 

J HI - sin 2 0 sm 2 l ) 


sin^ £7 sm* l 


sin i cos 0 sin 20 (c os i ~ 1 ) 

3 ( 4 sin2 0 sin2 1 “ X ) { X ” sin2 " cos i )} 2 (l - sin 2 0 sin 2 i) 


sin i cos 0 , „ * | sin 2 20(1 c os 1 ) 2 [ 

5 (4 sin 2 0 sin 2 1 - l) <1 - , . 2 Q . 2 «\ f sin 2 ^( cos 1 “ 1 ) 

° x ' 2 (1 " sin 2 0 sm 2 1) 


sin 20 (cos i - 1 ) 

+ cos 3 (\_ - K cos i sin i sin 0 (l - s in 2 0(1- cos i )} 2 T ““7 ” \ 

V 0 43 ' X V JS (l - sin 2 0 sin 2 i) 


cos i sin i sin 0 f sin 2 20(1 cosi) 2 | 

+ 9 .<1 - , " ~ sin 20 (cos i - 1 ) 

2 j 2 (l- sin 2 0 sin 2 i)J V ; 


+ sm x cos 0 (4 s . n2 Q s . n2 . _ jj |l - 2 (i _ s . ( n2 e sin2 ) i) } {l - s in 2 *(1 - cos i)} 

_ sin i^cos 8 ^ s . n2 e s . n2 . _ ^ _ s . n2 0(1 _ cos i)} ^ 0s . n2 .j l* (60) 
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The first two terms of the factor of sin 3(\ 0 - \ 43 ) in Equation 60 are odd functions about 8 = 77 . 

The third and fourth terms of the factor of sin 3(k 0 - \ 43 ) in Equation 60 are also odd functions 
about 6 - v, since they are controlled by sin 26 cos 8 - 2 cos 2 8 sin 8 which is "odd" about 8 = n. 

All the terms of the factor of cos 3 (a. 0 - a. 43 ) in Equation 60 are odd functions about 0 = W2 and 
3 V 2 . The first and second terms are controlled by sin 6 sin 26 = 2 sin 2 0 cos 6, and these are "odd" 
about 6 = 77/2 and 377/2. Thus, the orbit average of F t>43 (0) is zero and there can be no long-term 
contribution to the longitude drift of a 24-hour, inclined, circular- orbit satellite due to this gravity 
harmonic. 


Orbit-Averaged Effect of J 44 

The tangential force on the 24-hour inclined- orbit satellite at 6 from its ascending node, due 
to J 44 , is (see Appendix A) 

F T,44 (0) = F A,44 Sin a + F <6,44 COS a 

- k 44 -jsin 4 “ X- 44 ) sin a cos 3 <f> + cos 4 (A. ~ ^ 44 ) cos a cos 3 sin , (61) 

where 


M E 


k 44 = 420 (Ty J< 


Introducing the orbit geometry relations, Equations 29a through 29d causes Equation 61 to become 


F T,44 ^44 S sin4 ( A '0 ^“ 44 ) 


cos i (l - sin 2 9 sin 2 i) cos 4 AA. - sin 2 6 (l - sin 2 6 sin 2 i) sin 4A\J 

) cos 4AA.j^j- * (62) 


r * n 2 * 

+ cos 4(^ 0 ~ A 44 ) cos i (l - sin 2 9 sin 2 i) sin 4A\ + ~ sin 29 (l - sin 2 9 sin 2 i 


Using the longitude excursion relations. Equations 28 and 29, causes Equation 62 to become 


F x> 44 (0) = k 44 < S in4 (\ 0 -\ 44 ) 


cos i (l - sin 2 9 sin 2 i) 


2 [l - sin 2 (9(1 - cos i )] 2 [sin 2#(cos i - 1 )] 


(l - s in 2 6 s in 2 i ) 2 


. . f Tsin 29(cos i - 1 ) "1 f sin 2 20( 1 - cos i ) 2 

- sin 2 i sin 20 (l - sin 2 0 sin 2 i) {l - sin 2 6(1 - cos i )} ^ 7 — 7 - 77 — f ^ 1 - 


+ cos4(A. 0 -\, 4 ) 


sin 2 i sin 2 6 . _ ^ x 1 

+ ~ 2 (l - sin 2 9 sin 2 i) *s 1 - 


(l - sin 2 6 sin 2 i) J 1 2 (l - sin 2 6 sin 2 i) 


- f sin 20 ( cos i - 1 ) *) f sin 2 26 ( 1 - cos i ) 

i (l “ sin 2 9 sin 2 i) 2 {l - sin 2 6( 1 - cos i )) < : “ c > <1 - — : 

V *■ y/ \(l-sin2Ssin2 i)J \ 2(1 -sin 2 6 sin 2 i 

2 [l - sin 2 1 8(1 ~ cos i )J 2 [sin 26 (cos i - 1 )] | ^ (63) 

(l - sin 2 6 sin 2 i) 2 
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The factor of cos 4 (\ 0 - \ 44 ) in Equation 63 is an odd function about Q - v. Thus, it orbit averages 
to zero and has no relevance on the long term drift of the satellite. 

Considering only the relevant sin 4 (a. 0 - \ 44 ) term, Equation 63 becomes 


n'je) 


Ft , 44 (^1 " k 44 sin 4(A 0 - A 44 } -jcos i (l - sin 2 5 sin 2 i) - sin 2 i sin 29 [l - sin 2 9(1 - cos i)] sin 25(cos i - 1 )j- 

k 44 sin 4 (a.„ - A_ 44 ) r 

+ 2 (l - sin 2 9 sin 2 i) \ 4 COS 1 t 1 " sin2 ~ cos O] 2 [sin 25(cos i - 1 )] 2 

i 4 25(1 - cos i) 3 [l - sin 2 5(1 - cos i)]j- (63a) 


sirr i sm 


- k 44 sin4(^ 0 ~A 44 ^ ^A' (5) 

[l - sin 2 5( 1 - cos i )] [sin 25(cos i - 1 )] 2 [- 4 cos i [l - sin 2 5( 1 - cos i )} + sin 2 i sin 2 25(cos i - 1 )]1 /gov.\ 
+ 2 (1- sin 2 5 sin 2 i) 


~ k 44 sin 4(^ 0 -X 44 ) ‘SA' (5) 


l) 

2 [l - sin 2 5(1 - cos i)] [sin25(cos i - l)] 2 [- \+cos i jcos i +(l - sin 2 5sin 2 i)| - J+^osI sin 2 5(l - sin 2 5) Lg 3c ) 

[l - sin 2 6 sin 2 i} J 


B' ( 0 ) 


f *2 cos i [l - s in 2 6( 1 - cos i )] [s in 2<9(cos i ~ 1 )] 

k 44 sin4(A 0 -A 44 ) ^A' (5)- - ! + cos i 




2 [l ~ sin 2 6(1 ~ cos i )] [sin 2#(cos i - 1 )1 2 r 0 

/ 0 ~~~ 0 \ |l - sin 2 i (l - sin 2 i sin 2 6 ) - sin 4 i sin 4 6 f (63d) 

( 1 + cos 1 ) (1 - sm 2 6 sm 2 i) c v J 


= k 44 sin4(A 0 -X. 44 ) [A' (5) +B' (5)] 

k sin4ft fl -\ 4 ) rf A ^ ^ 

+ (1 + cosi) j^2 sin 2 i [l - sin 2 6( 1 - cos i)] [sin 2#(cos i - 1)] 2 


c' ( 6?) 



D' (6) 

= A___ 

sin 2 6( l - cos i )] [sin 2<9(cos i 

*v 


l)] 2 [l + sin 2 6 sin 2 i]^- 


(63e) 
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Thus, from Equation 63e 


F t>44 W) = k 44 sin 4 (^ 0 -X- 44 ) {A' < 0 )+B'( 0 )+C' ( 0 )+D' ( 0 )} 


= k 44 sin4(X 0 -\ 44 ) <Jcos i - cos l * l n 1 _ sin 2 i(cos i - l) ifj - sin 2 i( 1 - cos i) 2 (f 


i_1) (^ ) _ ! 


2 cos i( 1 - cos i ) (cos i - 1 ) 2 (^ ) 2 cos i(cos i - 1 ) 


(2) 


(1 + cos i ) 


( 1 + cos i ) 


2 sin 2 i(cos i 


( 1 + cos i ) 




2 sin 2 i(cos i-l) 2 (l~ cos i) l 4 


( 1 + cos i ) 


( cos i ” 1 ) 2 ( X 

+ hr 


(cos i + 1 ) \4 


2(1“ cos i ) (cos i - l) 2 2 sin 2 i(cos i “ l) 2 


(1 + cos i ) 


( 1 + cos i ) 


( 1 - cos i ) sin 2 i(cos i - l) 2 Ug 


(1 + cos i ) 


(64) 


, v . j cos 1 , _ . v sin‘ 1 , „ 

- k 44 sin4(^ 0 ~k 44 ) (2 - sin 2 i) + -5 (1 - cos 1) 


1 - 


( 1 - cos i ) 


cos i(cos i “ l) 2 

1 - cos i 

sin 2 i(cos i - l) 2 

( 1 “ cos i ) 

r 1 + cos i 

2 “ 1 

1 1 + cos i 

1 “ 2 


(cos i - l) 2 

( 1 ” cos i ) 

sin 2 i 

5 sin 2 i 

cos i + 1 

1 ~ 2 

+ 2 ■ 

16 ( 1 " cos 1 )/ 


(64a) 


. . I cos 1 / _ v s in 1 

k 44 sin4(\ 0 ~\ 4 ) 2 — (l +cos 2 i) + 4 (1 - cos i) (1 + cos i) 


(cos i “ 1 ) 2 
1 + cos i 


(1+cosi) + 2 (1 + cosi) ~ 


( 1 + cos i ) 


1 2 4 5 sin 2 i(l - cos i) 


16 


(64b) 
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= k 44 sin4K-X 4 J ( 1 + cos 2 i)+ (l ~ C ° S 4 2i) - 


(cos i - l) 2 


. * (cos i - l) 2 (1 - cos i) 

(1 + cos i) + (l - cos 2 i) - — - (3 + 5cosi) 


(64c) 


, . /. x \ J cos i cos 3 i 1 cos 2 i cos 4 i cos 3 i 0 . 

k 44 sin 4 [K 0 -\ 44 ) i — 2 + — T~ + T “ ~~2 — + ~4 2 + cos 1 


cos 1 COS 7 1 


+ COS' 5 1 


,2 j (3+5 cos i ) 


16 


( 1 - 3 cos i + 3 cos 2 i - cos 3 i)|^ (64d) 


k 44 sin4 


( x o ^ 44 ) {ie 


4 cos i 6 cos 2 i , 4 cos* 
16 T 16 + 16 + 16 


i , cos 4 i L 

ie ; 


(64e) 


- k 44 sin4(^ 0 ~^ 44 ) 


( 1 + cos i ) 4 

16 


(64f) 


But k 44 sin 4(X- 0 -X 44 ) = F 44 sin4(^ 0 -^ 42 ) , is the tangential perturbation force on the 24-hour 
equatorial satellite at longitude \ 0 . Thus, the orbit-averaged perturbing force on a 24-hour, 
inclined circular orbit satellite whose ascending node is at is just F(i) 44 times the perturbing 
force on the equatorial 24-hour satellite at ^ 0 , due to J 44 , where 

( 1 + cos i ) 4 

F(i) 44 = 16 • ( 65 ) 


Summary 

The orbit- averaged, long-term drift of the inclined circular orbit, 24-hour satellite in an 
earth gravity field through fourth order is derivable from 


A. = 


■ (i2tt 2 ) 


q n 

LZ> nm sinm(k -k nn ) F(i) nm , rad/(sid. day) 2 


n = 2 nv= 1 
For n-m, even 


where the F nm are given in Equation 2, and 


( 66 ) 


( 1 + cos i ) 2 
F(i) 22 = 4 


(67a) 
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(67b) 


and 


F(i) 


31 


(1+cosi) 5 sin 2 i(l + 3 cos i ) 
2 8 


F(i) 
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( 1 + cos i ) 3 
8 


w N _ (1+cosi) 2 7 sin 2 i cos i( 1 + cos i ) 

F( 1 ) 4 2 - ” A a 


F(i) 


44 


( 1 + cos i ) 4 
16 


( 47r a 0 y -1 . / . v length units 

g s sid. day 


n = 2 m= 1 
For n~m, even 


FIRST INTEGRAL OF THE LONG-TERM LONGITUDE DRIFT 


(67c) 

(67d) 

(67e) 


( 68 ) 


The "trajectory solution" of \ and a as a function of time for the 24-hour satellite drift in the 
fourth-order gravity field, stemming from Equations 66 and 68, cannot be expressed in a closed 
form with elementary functions. However a first "energy integral" of Equation 66 can be found 
which will prove valuable in the applications. 

The common separation technique used in Reference 5 yields: 


•• = dA. d (A.) 2 _ d(\) 2 

dt “ 2A. dt = 2d\ * 

Substitution of Equation 69 into Equation 66 enables one to separate variables giving 


(69) 


d(M 2 = — if^2 ZZ] F "« sinm ^~ X "-) F(i) ™ dX • ( 7 °) 

Integration of Equation 70 under the assumption that the only variables are k, the geographic longi- 
tude of the ascending node of the satellite, and k, its drift rate, gives 

(k) 2 = C 0 + cos [rad/s id . day] 2 . (71) 
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In the applications, C 0 in Equation 71 can be considered to be the mean long-term squared drift 
rate of the 24- hour satellite in the drift period under consideration. C 0 is evidently a function of 
the initial elements of the orbit of the satellite as well as all the earth gravity constants and the 
initial positions of the sun and moon (see the section, COMPARISON OF THE THEORY WITH 
NUMERICALLY INTEGRATED TRAJECTORIES). 

SIMPLE CLOSED-FORM SOLUTION TO THE LONGITUDE DRIFT OF A 24-HOUR 
SATELLITE FOR PERIODS CLOSE TO SYNCHRONOUS 

If the period of a 24-hour satellite is sufficiently close to synchronous, Equation 66 describes 
the motion of the mean geographic longitude of this satellite of small orbit eccentricity, due to 
earth gravity through fourth order (see DISCUSSION). 

With the definitions of the constants A nm given in Equations 76, 85, 94, 102, and 105, Equation 
66 becomes 


4 n 

= 22 22 A " m sinm ^ ■ ( 7ia ^ 
n = 2 m - 1 
n - m, even 

If the problem is restricted to small excursions in mean longitude from an initial 24-hour satel- 
lite configuration, Equation 71a may be integrated twice to give the long-term drift explicitly in 
a simple closed form. 

Let this small excursion AA. be given from an initial configuration at mean longitude (or longi- 
tude of the ascending equator crossing) >^ 0 , so that 

A. = A.© + AA. . (71b) 

Providing m/Sk is sufficiently small, when Equation 71b is substituted into Equation 71a and it is 
expanded, the small excursion, long-term acceleration due to fourth- order earth gravity becomes 

^ = T.T Sinm (\ 0 -X. n J A n " C ° S m K _ \>l m) • ( 71C ) 


With Equation 71c rearranged, the differential equation of long-term longitude drift for small 
excursions is 


A^ + A^ 


-m. 


mcosm(\ 0 -\ niii ) 


V -1 V 


/ ,7 , A nn, sinm ( X Q-\,J 


(71d) 
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A trial solution to the homogeneous equation of Equation 7 Id is 


AX h - Csinoit + D cos cot 


(71e) 


Equation 71e satisfies the homogeneous equation of Equation 71d only if 


CO - + 


^~!7~^ A nn. m COS * ( X 0 " A nm) 


-1 1/2 


(7 If) 


The particular solution to Equation 7 Id is 


AX 




XiI] A nm mCOS m 


(71g) 


The complete solution to Equation 7 Id is thus 


AK = A\, + A\ 

h p 


C sinwt + D cos cot 


Jj^j A nm Sinm (\ 0 ~Km) 

££ A nm mCOSn 1 ( X 0 ~\ lm ) 


(71h) 


An initial value problem corresponding to Equation 71d may be specified by giving the initial drift 
rate AX 0 (small) at time t = 0, when AX = o . 

Thus, let 


AX - 0 , when t - 0 , 


and 


AX - AX qj when t = 0. 


This specification in the complete solution, Equation 71h, determines the coefficients c and D: 

D = 

ZL A n m m c Osm(\ 0 - A nm) 
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and 


C 


AX n 


Cl) 


The small excursion from an initial near 24-hour configuration at \ 0 and t = 0, is thus given from 
Equation 71h as 


AX 


^0 

s in cot + 


1] ZXm Sinm K ~0 

2] 2] A n"> mCOSm ( A -0"^n m ) 


[cos Olt ~ l] 


(7 li) 


It may now be asked under what conditions will the mean longitude excursion of a 24-hour satellite 
be negligibly small for all time under the influence of earth gravity through fourth order. Such 
potentially stable mean longitude positions for the 24- hour satellite are to be found at those 
longitudes X 0 , where (from Equation 71i) 

4 n 

A n. sini "K-0 = 0 ' ( 71 j) 

n = 2 m= 1 
n~m, even 

For the equatorial 24-hour satellite, the general stability criterion for 24-hour satellites, Equa- 
tion 71j, was first stated in Reference 5. There it was shown that for the M real earth” there are 
four such points around the equator, where long term stability (on the order of years) can be 
achieved for this satellite in the presence of earth gravity through fourth order. Of these four 
longitudes, two are dynamically stable and two are only statically so. The first two points fall at 
the two relative minima of the longitude drift potential for the 24- hour equatorial satellite. The 
two longitudes only statically stable fall at the two relative maxima of the longitude drift potential 
(see Figure 5). In Reference 6 it was shown that all of these stability points, in fact the entire 
drift regime for the very high 24-hour equatorial satellite of the M real earth,” are strongly dom- 
inated by the second order tesseral harmonic J 22 that is associated with the ellipticity of the 
equator of the earth. A comparison of the relevant inclination factors F(i) nm , Equations 67, will 
show that equatorial ellipticity dominates the drift regime due to the earth of the inclined 24-hour 
satellite to an even greater degree for all inclinations less than 90°. However, for the inclined 
satellite, the total strength (for example the absolute difference between its relative maxima and 
minima) of this potential declines roughly in proportion to (1 + cos i) 2 / 4. Thus, for highly inclined 
24-hour satellites, sun and moon perturbations may have appreciable influence on long-term 24- 
hour satellite drift. 

Returning to Equation 71i, it may be seen that the criterion for the type of long term stability 
for satellite placements near \ 0 with zero initial drifts, where sinm (\> - \ J io > is 

established by the character of the ''frequency" term at the placement longitude \ 0 . If 
-EE A nm m cos m (x„ - K n Jj at the "stable" longitudes A 0 (from Equation 71j)is greater than zero, 
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then (from Equation 7 If) is a real number and the character of the long term drift in the neigh- 
borhood of k 0 is harmonic or self- limiting. These longitudes are the dynamically stable ones, 
and they have been shown (References 1 and 5) to be at or near the longitudes of the minor equa- 
torial axes for the "real earth" (\ 22 ± 90°) . On the other hand, if - A nm m cos m (^o " ^ nm ) is less 
than zero at or near a "stable" longitude k Q9 co is an imaginary number and the character of the 
long term drift in the neighborhood of k Q will be exponentially divergent. These longitudes are the 
merely statically stable ones and they have been shown (References 1 and 5) to be at or near the 
longitudes of the major equatorial axis for the "real earth" ( \ 22 ) . 

As to the excursion range of validity of Equation 71i, numerical evaluation of the exact elliptic 
integral of J 22 -drift (Reference 6, Appendix D), in comparison with Equation 7 li indicates it will 
predict excursions to the order of mA^ = 20°, with an order of accuracy of about 1%. Since J 22 
dominates the drift of the 24-hour satellite at almost all longitudes, this implies that Equation 71i 
is generally useful for predicting long term excursions of the 24- hour satellite from an initial 
mean longitude \ 0 up to about 10 degrees, providing A\ 0 is sufficiently small (see DISCUSSION on 
this latter point), and providing sun and moon perturbations are negligible. 

To illustrate the accuracy of Equation 71i for predicting long term high order drift in the 
neighborhood of the major equatorial axis of the earth where the influence of j 22 is small, a com- 
parison of Equation 7 li with a numerically integrated third order earth, sun, and moon gravity 
drift is presented in Table 4, Case 2. The earth gravity field used is a recent weighted average 
estimate due to W. M. Kaula (Appendix A). 

The coefficients and initial orbit parameters in Table 4, Case 2, give the following: 

A 22 = 20.7949 X 10“ 6 , rad/sid. day 2 , 

A 3I = -.254774 x 10“ 6 , rad/sid. day 2 

A 33 = 2. 14239 x 10“ 6 , . rad/sid. day 2 , (71k) 


LE 


m cos m 

n m 


i k o~K m ) = 35.5169 X 10' 6 , rad/sid. day 2 


(711) 


E L Anm sinm ^ 


-6.4100 x 10 6 rad/sid. day 2 


(71m) 


Equation 71 1 substituted in Equation 7 If gives 


a> = i (5 . 95961 x 10' 3 rad/sid. day) 
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Thus Equation 7 11 becomes [noting cos iat = coshat, sin iat = i s inhat] 


A\ 


. 02289°/solar day sinh [s. 95961 x 10 3 t] 
o v t 

5.95961 x 10“ 3 rad/sid. day 


. 180477 


jcosh [5. 95961 X 10“ 3 t] 



= (.0668525) sinh [5. 97593 x 10" 3 t] - .180477 jcosh [5. 97593 x 10' 3 t] - lj radians , (71n) 

with t (in the final form of Equation 7 In) in units of solar days. Equation 71h is evaluated in the 
last column of Table 4. The discrepancy between the theoretical and the numerically integrated 
drift illustrates the relatively large influence solar-lunar gravity can have on the drift of a 24- 
hour satellite near the longitude of the axes of the elliptical equator of the earth (see Figure 5b). 


COMPARISON OF THE THEORY WITH NUMERICALLY INTEGRATED TRAJECTORIES 

Tables 1, 2, 3 and 4 and Figures 2, 3, 4 and 5 give data from numerically integrated 24-hour 
particle trajectories about the earth in various perturbing gravity fields. 

The trajectories were computed utilizing an earth gravity field as represented in Appendix A. 
The Gaussian gravity constant of the earth used for all the trajectories was 


= 3.9860319 x 10 s km 3 /sec 2 . 

The mean equatorial radius of the earth used for all the trajectories was 


R 0 = 6378.165km 


Tables 1, 2 and 3 give trajectories in earth gravity perturbing fields due to the tesseral harmonics 
j 22 , J 31 and J 33 respectively, acting alone on a 24-hour satellite of zero and 60 degree orbit 
inclination. 


The J nm longitude harmonic constants were chosen larger than they are in reality for these 
cases to emphasize the effects of the perturbations in a reasonable period of time. 

The additional earth gravity constants used in these trajectories were (Tables 1, 2, and 3; il- 


lustrated in Figures 2, 3 and 4) 


J 2 2 

-6.0 X 10' 6 , 

J 3 1 

-100.0 x 10~ 6 

J 33 

-10.0 X 10' 6 


^22 = 

-21.0° , 




II 

r*> 

-156.0° , 


> 

(72) 

^33 = 

-36.0° . 

y 
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Table 1 


24-hour Satellite Drift (J 22 perturbation only). 
(j„ = -6.0 x10-6,a = -21.0°) 


INJECTION CON DITIONS: CASE 1 
Inclination: 0.0° 

Eccentricity: .149 x 10' 7 

Initial Drift Rate: -,00566°/day at -66.00567° 


INJECTION (CONDITIO NS: CASE 2 
Inclination: 60.0° 

Eccentricity: .149 x 10“ 7 

Initial Drift Rate: .00322°/day at -66.00323° 


Time 

Numerically 

Calculated 

Numerically 

Calculated 

Theoretical 

Time 

Numerically 

Calculated 

Numerically 

Calculated 

Theoretical 

(days from 
injection) 

Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Longitude 
Drift Rate 
(°/day) 

Drift Rate 
(°/day) 

(days from 
injection) 

Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Longitude 
Drift Rate 
(°/day) 

Drift Rate 
C/day) 

0.00000 

-66.00001 



0.00000 

-66.00001 



(1.00000) 

(-66.00567) 

-.00566 

-.00566 


(-66.00323) 

-.00322 

-.00322 

2.00000 

-66.01132 



1.99454 

-66.00644 



13.9687 

-66.54857 



13.96261 

-66.30873 



(14.9635) 

(-66.63223) 

-.08409 



(-66.35585) 

-.04724 


15.9583 

-66.71588 



15.95742 

-66.40297 



29.9271 

-68.51610 



29.92199 

-67.41438 



(30.9271) 

(-68.68971) 

-.17361 



(-67.51189) 

-.09775 


31.9271 1 

-68.86331 



31.91706 

-67.60939 



43.8958 

-71.40683 



43.88826 

-69.03966 



(44.8958) 

(-71.65802) 

-.25119 



(-69.18039) 

-.14106 


45.8958 

-71.90921 



45.88360 

-69.32112 


i 

57.8750 

-75.37480 



57.85623 

-71.27551 



(58.8698 

(-75.69993) 

-.32682 

-.32701 


(-71.45938) 

-.18428 

-.18520 

59.8646 

-76.02505 



59.85180 

-71.64325 


i _ 


Table 2 


24-hour Satellite Drift (J 31 perturbation only). 
(j 31 = -100.0x10-6,^3! = -156°) 


INJECTION CON DITIONS: CASE 1 
Inclination: 0.0° 

Eccentricity: .149 x 10" 7 

Initial Drift Rate: -.00358°/day, at -66.00359° 


INJECTION CONDITIONS: CASE 2 
Inclination: 60.0° 

Eccentricity: .149 x 10" 7 

Initial Drift Rate: .00 146° /day at -65.99855° 


Time 

Numerically 
Calculated 
Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Numerically 

Calculated 

Theoretical 

Time 

Numerically 
Calculated 
Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Numerically 

Calculated 

Theoretical 

(days from 
injection) 

Longitude 
Drift Rate 
(°/day) 

Drift Rate 
(°/day) 

(days from 
injection) 

Longitude 
Drift Rate 
(°/day) 

Drift Rate 
(°/day) 

0.00000 

-66.00001 



0.00000 

-66.00001 




(-66.00359) 

-.00358 

-.00358 


(-65.99855) 

.00146 

.00146 

2.00000 

-66.00716 



1.99451 

-65.99709 



13.96875 

-66.34588 



13.96135 

-65.85450 




(-66.39862) 

-.05302 



(-65.83231) 

.02226 


15.95833 

-66.45136 



15.95578 

-65.81011 



29.92708 

-67.58556 



29.91624 

-65.33211 




(-67.69498) 

-.10942 



(-65.28606) 

.04619 


31.92708 

-67.80440 



31.91051 

-65.24000 



43.89583 

-69.40947 



43.87586 

-64.56675 




(-69.56747) 

-.15883 



(-64.50049) 

.06645 


45.88542 

-69.72547 



45.87005 

-64.43423 



57.86458 

-71.92084 



57.83472 

-63.51238 




(-72.12761) 

-.20785 

-.20814 


(-63.42540) 

.08724 

.08769 

59.85417 

-72.33437 



59.82879 

-63.33842 
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Table 3 


24-hour Satellite Drift (J 33 perturbation only). 
(J33 = -10.0 x 10' 6 , ^33 = -36.0°) 


INJECTION CONDITIONS: CASE 1 | 

INJECTION CONDITIONS: CASE 2 


Inclination: 0.0° 



Inclination: 

60.0° 



Eccentricity: .149 x 10 7 


Eccentricity: .149 x 10 7 



Initial Drift Rate: - 

.01066° /day, at -66.01067° 

Initial Drift Rate: -.00454° /day, at -66.00454° 


Numerically 

Calculated 

Numerically 



Numerically 

Calculated 

Numerically 


Time 

Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Calculated 

Theoretical 

Time 

Longitude of 
the Ascending 
Equator Crossing 
(degrees) 

Calculated 

Theoretical 

(days from 
inj ection) 

Longitude 
Drift Rate 
(°/day) 

Drift Rate 
(Vday) 

(days from 
injection) 

Longitude 
Drift Rate 
C/day) 

Drift Rate 
(°/day) 

0.00000 

-66.00001 



0.00000 

-66.00001 




(-66.01067) 

-.01066 

-.01066 


(-66.00454) 

-.00454 

-.00454 

2.00000 

-66.02133 



1.99455 

-66.00906 



13.99687 

-67.03684 



13.96297 

-66.43768 




(-67.19582) 

-.15898 



(-66.50455) 

-.06704 


15.99687 

-67.35480 



15.95790 

-66.57142 



29.93750 

-70.75079 



29.92364 

-68.00616 




(-71.07700) 

-.32621 



(-.68.14449) 

-.13865 


31.93750 

-71.40321 



31.91893 

-68.28281 



43.91667 

-76.14489 



43.89178 

-70.30889 




(-76.60462) 

-.46213 



(-70.50875) 

-.20030 


45.90625 

-77.06434 



45.88742 

-70.70861 



57.89583 

-83.30173 



57.86230 

-73.46010 




(-83.87285) 

-.57112 

-.57160 


(-73.71939) 

-.25981 

-.25953 

59.89583 

| -84.44396 



59.85826 

-73.97868 




( J 22 ~ 


Time 
(Solar 
days from 
injection) 


1.0 

24.9 

50.85 

74.8 

100.75 

124.7 

148.6 


Table 4 


24-hour Satellite Drift in a Third Order Earth, Sun and Moon Gravity Field. 

-1.51 x 10 -6 , A -22 = -15. 5°, J 31 = -1.51 x KT 6 , ^ 31 = 0.0°, J 33 = -.149 x 10' 6 , X 33 = 22.8°) 


INITIAL CONDITIONS: CASE 1 


INITIAL CONDITIONS: CASE 2 


Inclination: 32.8° 



Inclination: 

32.8° 


Eccentricity: .45 x 10 7 , (a s /R 0 ) 

= 6.6107211) 

Eccentricity: .745 x 10" 8 , (a s /R 0 ) = 6.6107208 

Longitude: -66.00001° 


Longitude : - 

-22.47719° 


Longitude Rate: 

00554° /day 


Longitude Rate: .02289° /day 


Numerically 
Calculated 
Longitude of 
the Ascending 
Equator Crossing 

Numerically 
Calculated 
Longitude 
Rate (°/day) 

Theoretical 

Longitude 

Rate 

(Vday) 4. 

Time 
(Solar 
days from 
injection) 

Numerically 
Calculated 
Longitude of 
the Ascending 
Equator Crossing 

Theoretical 

Longitude 

(degrees) 

(degrees) 




(degrees) 


-65.99449 

.00554 

.00544 

0.00000 

-22.47719 

-22.47719 

-66.12636 

-.01827 


22.93488 

-22.00060 

-22.048 

-66.91847 

-.04323 


48.86178 

-21.67153 

-21.79 

-68.24578 

-.066799 


74.78919 

-21.55145 

-21.76 

-70.34043 

-.093496 

-.09415 

98.72270 

-21.62043 

-21.94 

-72.86641 

-.11874 


122.65685 

-21.89455 

-22.32 

-75.94402 

-.13992 

-.13970 

148.58637 

-22.47249 

-22.97 
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■Ill III I I 


LONGITUDE EAST OF GREENWICH 


-76.0° 


-75.0° 


-74.0° 

“I 


THE LONGITUDE AND DRIFT RATE REFERS TO THE LONGITUDE 
OF THE ASCENDING EQUATOR CROSSING OF THE TRAJECTORY 


NOTES: 

EARTH LONGITUDE GRAV- 
ITY CONSTANTS USED: 

J 22 = -6.0 x 10~ 6 
X 22 - -21.0° 

(SEE APPENDIX A) 


-Kl 


58.87 



i = 0°, ACTUAL DRIFT VS. DRIFT RATE- 

THEORETICAL DRIFT VS. DRIFT RATE 


DAYS FROM INJECTION 
ACTUAL DRIFT , 

NUMERICALLY CALCULATED 
BY AN ENCKE INTEGRA- 
TION METHOD. INITIAL 
ECCENTRICITY OF THE 
24 HOUR ORBITS = .149x 10' 7 . 

THEORETICAL DRIFT . 

FOR A CIRCULAR 24 HOUR 
ORBIT. 


H -.2 


H-.3 


Figure 2—24-hour satellite drift due to J 22 only. 


THEORETICAL DRIFT VS. DRIFT RATE 


LONGITUDE WEST OF GREENWICH 


73.0° 


72.0° 

_L 


71.0' 

_L 


70.0° 

X 


69.0° 


68 . 0 ° 


58.86 



l = 0°, ACTUAL DRIFT VS. DRIFT RATE- 
44. 87^ 


THEORETICAL DRIFT VS. DRIFT RATE 



58.83 


+ 0.1 


.89 

i = 60. 0°, ACTUAL 
DRIFT VS. 

DRIFT RATE . 


0 . 0 ° 


NOTES: 


EARTH LONGITUDE 
GRAVITY CONSTANTS USED: 

J 31 = -100.0 x 10~ 6 
X 31 = -156.0° 

(SEE APPENDIX A) 

DAYS FROM INJECTION 

50 3 -ACTUAL DRIFT 

NUMERICALLY CALCULATED 
BY AN ENCKE INTEGRATION 
METHOD. INITIAL ECCENTRI- 
CITY OF THE 24-HOUR ORBITS 
* . 149 x 10~ 7 

.^ THEORETICAL DRIFT . 

FOR A CIRCULAR 24 -HOUR ORBIT. 
THE LONGITUDE AND DRIFT RATE _| 
REFERS TO THE LONGITUDE OF 
THE ASCENDING EQUATOR 
CROSSINGS OF THE TRAJECTORY. 


- 0.1 


- 0.2 


Figure 3—24-hour satellite drift due to J 31 only. 
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DRIFT RATE (degrees /day) DR1FT RATE (degrees / day) 



-84.0° 


-82.0° 

I 


-80.0° 

I 


LONGITUDE EAST OF GREENWICH 

-78.0° -76.0° -74.0° -72.0° 

I 


-70.0° 


- 68 . 0 ° 


- 4 - 0.1 


H -0.2 



H -0.3 


NOTES: 

EARTH LONGITUDE GRAVITY 
CONSTANTS USED: 

J 33 = - 10.0 x 10 -6 
A 33 = -36.0° 

(SEE APPENDIX A) 

DAYS FROM INJECTION 


H-0.4 


58.90 


THEORETICAL DRIFT VS. DRIFT RATE 


ACTUAL DRIFT , NUMERICALLY 
CALCULATED BY AN ENCKE 
INTEGRATION METHOD, 

INITIAL ECCENTRICITY OF 
THE 24-HOUR ORBITS = .149x 10" 

THEORETICAL DRIFT FOR 
A CIRCULAR 24-HOUR ORBIT. 

THE LONGITUDE AND DRIFT RATE 
REFERS TO THE LONGITUDE OF 
THE ASCENDING EQUATOR 
CROSSING OF THE TRAJECTORY 


H -0.5 


H -0.6 


Figure 4 — 24-hour satellite drift due to J 33 only. 


The "synchronous " semi-major axis a s , assumed in the theory to be constant for all the orbits of 
these trajectories, is taken as the semi- major axis at "injection" of the 24- hour satellite, and is 


a s = 42,164.27 km . 


Thus, in the theory, the constant (R 0 / a s ) is 


(VO = (fiiM 1 #) = 1512694 • 


(73) 


the constant (R 0 / a s ) 2 is 


(R 0 /a s ) 2 = .02288243 . 


(74) 
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DRIFT RATE (degrees /day ) 



LONGITUDE OF THE ASCENDING EQUATOR CROSSING - WEST OF GREENWICH 


76 ° 


73 

I 


LONGITUDE WEST OF GREENWICH 
72° 71° 70° 

I 



THE LONGITUDE AND DRIFT RATE REFERS TO THE LONGITUDE 
OF THE ASCENDING EQUATOR CROSSING OF THE TRAJECTORY. 

NOTES: 

EARTH GRAVITY CONSTANTS USED: 

J 20 = 1082.30 x 10~ 6 , J 30 = -2.30 x 10“ 6 , 

J 40 - -1.80 x 10“ 6 , J 22 = -1-51 x 10" 6 , 

J 31 = -1.51 x 10 -6 , J 33 = - .149 x 10 " 6 
X 22 = ~ 15.5° , X 3l = 0.0°, X 33 = 22.8° 

(SEE APPENDIX A) 


148.6 


THEORETICAL DRIFT VS. DRIFT RATE 


DAYS FROM 
INJECTION 
(30.0 JAN. 1965) 


100.75 _ 

THEORETICAL DRIFT ; 

S FOR A CIRCULAR 24 -HOUR 
ORBIT OF INCLINATION 

ACTUAL DRIFT . 1 = 32 - 8 °' ,N AN EARTH FIELD 

NUMERICALLY CALCUALTED ONLY. 

50.0 BY AN ENCKE INTEGRATION- 
METHOD. INITIAL ORBIT 
ECCENTRICITY: .45 x 10 


H “0.1 


Figure 5a— 24-hour satellite drift in a third order earth, sun and moon gravity field. 


23° 


ACTUAL DRIFT ; CALCULATED 
BY AN ENCKE INTEGRATION 
METHOD. INITIAL ORBIT 
ECCENTRICITY: .745 x I0‘8 


THEORETICAL DRIFT ; 

FOR A CIRCULAR 24-HOUR 
ORBIT OF INCLINATION 
i = 32. 8°, IN AN EARTH 
FIELD ONLY. 


/ 

/ 

/ 

/ 

/ 

/ 

/ 


22 ° 


21 ° 




TIME FROM INJECTION, 31.0 JAN. 1965 (solar days) 


Figure 5b— 24-hour satellite drift in a third order earth, sun and moon gravity field. 
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DRIFT RATE (degrees /day ) 



and the constant (R 0 /a s ) 3 is 


( R o/ a s) 3 = 003461411 . (75) 

The Theory for 24-Hour Satellite Drift in a "J 2 2 -only” Perturbation Field 

When Equations 71, 2 and 67a are combined and it is noted from Equation 8 that g s = M E /( a s ) 2 , 
Equation 71 becomes for ”J 22 -only" drift 

(M 2 = C 0 - A 2 2 cos 2 (A. - 2 ) , (76) 


where 


(cos i + l ) 2 /r7r7 v 

A 22 = ~ 72?r2 ( R o/ a s) 2 J22 4 (rad/sid. day ) 2 . ('*) 

(Compare with Reference 6, Equation 36.) 

The corresponding longitude acceleration during ’’ j 22 -only" drift is given from Equation 65 as 

A. = A 22 sin 2 ~ ^ 22 ) . (77a) 

The coefficient A 22 in Equation 77a is always positive since all the tesseral j nm ’s in the gravity 
potential form of Appendix A are arbitrarily assigned as negative numbers. Equation 77a es- 
tablishes a one-dimensional potential regime for the long term movement of the mean longitude 
of the circular orbit 24-hour satellite. This potential has a shape derivable from the right side of 
Equation 77a which is the effective ’’force” moving the longitude of the satellite. For a one- 
dimensional potential, the force is given as the first derivative of the potential. The one- 
dimensional (longitude) potential energy of the satellite (the negative of the potential) associated 
with the longitude perturbation expressed by Equation 77a is thus (assigning the potential constant 
as zero) 


(P ' E,) *22 ” + ~^ COs 2 (^— \ 2 2 ) 
(see Figure 6) 


(77b) 


The longitude potential energy of the satellite s from Equation 77b thus has an absolute and a 
relative maximum at K = X 22 , and at k = ^ 22 + 180°. These are evidently points of unstable 
equilibrium where the longitude perturbation force from Equation 77a is momentarily zero. The 
potential energy due to J 22 and A 22 (from Equation 77b) has an absolute and a relative minimum 
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GREENWICH GREEN - 



Figure 6— 24-hour safe! life Jongifude-driff pofenfial 
energy and drift regimes due to J 22 . 


at k -k 22 = ± 90 ° , or k = ±90° + ^ 22 , which are at 
longitudes lying on the extensions of the equa- 
torial minor axis. These are evidently the only 
points of stable equilibrium for the placement 
of 24- hour satellites considering only the j 22 
perturbation. The "pendulum” dynamics as- 
sociated with Equation 77a has been discussed 
previously in Reference 6 in greater detail. 

If the initial longitude drift rate of the sat- 
ellite is A- 0 , with the ascending equator crossing 
at k 0 , then c 0 in Equation 76 is evaluated from 


(^o ) 2 C o A 22 cos 2 (X 0 ^22) 

giving 


C 0 - K) 2 + A 2 2 COS 2 (\>- X 2 2 ) 


(78) 


Substitution of Equation 78 into Equation 76 gives the theoretical drift rate of the 24- hour satellite 
at any subsequent longitude k 9 as 


k = ± |(\ 0 ) 2 + A 22 [cos 2(^o ~ k 22 ) ~ cos 2(k - ^ 22 )]| 1/2 rad/sid. day . (79) 

(The sign of the root in Equation 79 is to be taken according to the sense of the drift with respect 
to the longitude where k = 0 for the libratory regime (Reference 6).) 

From Equations 72, 74 and 77, A 22 is 


A 22 (i = 0°) = 97.56311 x 10 6 (rad/sid. day) 2 


(80) 


and 




A 22 (i=60°) = 97.56311 x 10" 6 x - — 4 


= 54.87925 x 10' 6 (rad/sid. day) 2 


(81) 
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For the case of the equatorial 24- hour satellite in Table 1, (Case 1) 


= -66.00567 , 


(^o ) 2 


00566°/ s olar day 
57. 29578°/ rad 


x .9972696 


so lar day 
sidereal day 


= 97.05385 x 10“ 10 (rad/sid. day) 2 


(82) 


Substitution of Equations 72, 80 and 82 into Equation 79 gives the theoretical drift rate for the 
equatorial 24-hour satellite in Table 1, Case 1 as 

k = - {97.054 x 10' 10 + 97.563 x 10" 6 [- . 00020 - cos 2 K + 21 . 0°)]j 1/2 (rad/sid. day). (83) 

See Table 1, Case 1 for an evaluation of Equation 83 at the respective ascending equator- crossing 
longitudes listed for the numerically calculated drift (see also Figure 2). 

For the case of the 60.0° inclined- orbit 24-hour satellite in Table 1, the theoretical drift 
formula corresponding to Equation 83 is 

k = - {31.412 x 10" 10 +54.879 x 10“ 6 [-. 00011 - cos 2(k + 21 . 0°)]| 1/2 (rad/sid. day) (84) 

(see Table 1, Case 2 and Figure 2). These cases make it clear that the theory for the drift of a 
24-hour circular orbit satellite due to J 22 , summarized in Equations 76 or 79, is essentially an 
accurate one. In the "Discussion" the probable limits of accuracy of the theory with respect to ec- 
centricity and period are stated more precisely. 


The Theory for 24-Hour Satellite Drift in a f, J 31 -only" Perturbation Field 

Combining Equations 71, 2 and 67b as in the first part of this section causes Equation 71 to 
become for "j 31 -only" drift 


M 2 - C 0 ” 2A 3 1 COS { k ~ k 3l) 


(85) 


where 


18" 2 (V a s) 3 J 


■■F 


5 sin 2 i s 
8 


*4 


(l + 3cos i ] [ (rad/sid. day) 2 . (86) 
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The corresponding longitude acceleration during ”J 31 -only” drift, is given from Equation 66 as 


k = A 0 


( k ~ k 3 1) 


(86a) 


Equation 86a establishes a one-dimensional potential energy for the longitude excursions of the 
24-hour satellite with respect to the J 31 perturbation. As before, this energy function is given as 
the negative of the first integral of the right hand side of Equation 86a (ignoring the integration 
constant): 


(P-E.)a 31 A 31 cos (k X 31 ) . 

(see Figures 7 and 8). 


(86b) 



Figure 7—24-hour satellite longitude-drift potential 
energy and drift regimes due to J 31 (A 31 < O, -41° < i 
< 41 |i| > 95.1°). 



Figure 8—24-hour satellite longitude-drift potential 
energy and drift due to J 31 (A 31 > O, 41° < li [ < 95. 1°). 


Depending on the inclination of the 24- hour satellite, two oppositely directed libratory longitude 
drift regimes are possible due to J 31 (see "NONRESONANT” INCLINATIONS FOR THE EARTH 
GRAVITY HARMONICS j 31 AND j 42 ). In the first, for |i| <41° and i between ±180° and ±95.1°, 
stable long-term librations can occur about k = \ 31 . In the second regime, for i between ± 41 ° and 
±95.1°, stable long-term librations can occur about k = k 3i +180° (see Figures 7 and 8 and Ap- 
pendix A). At the critical inclinations i = ±41°, +95.1° and ±180°, of course, there is no long-term 
drift perturbation at any longitude due to j 3I since A 31 = 0 for these inclinations. 

If the initial longitude drift rate of the satellite is k 0 with the ascending equator crossing at 
k 0 , then C 0 in Equation 85 is evaluated as in Equation 78 to yield, for ”j 31 -only” drift, 

(k) = ±|(\ 0 ) 2 +2A 31 |cos(\ 0 -A 31 ) -cos(\ '^ 31 )]| 1/2 (rad/sid. day) (87) 

(The sign of the root in Equation 87 is established from the sense of the drift with respect to the 
longitude where k = 0 for the libratory regime.) 
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From Equations 72, 75 and 85, 2A 31 is 


2A 31 (i = o°) = -122.9859 x 10' 6 (rad/sid. day) 2 , 


( 88 ) 


and 


2A 3 , (i = 60°) = 51.88468 x 1CT 6 (rad/sid. day) 2 . 

For the case of the equatorial 24- hour satellite in Table 2, (Case 1) 

\ 0 = -66.00359° 


and 


(89) 




(.00358) x (.9972696)1 2 
57.29578 


— = 38.82808 x 10" 10 (rad/sid. day) 2 


(90) 


Substitution of Equations 72, 88 and 90 into Equation 87 gives the theoretical drift rate for the 
equatorial 24-hour satellite in Table 2, Case 1 as 


K = - |38.828 x 10“ 10 - 122.99 x 10~ 6 [. 00006 - cos (\ + 156 . 0° )]j 1/2 ( rad/sid. day) . (91) 

See Table 2, Case 1 and Figure 3 for an evaluation of Equation 91 at the maximum longitude drift 
excursion of the satellite in the numerically computed trajectory. 

For the case of the 60.0° inclined- orbit, 24-hour satellite in Table 2, the theoretical drift 
formula corresponding to Equation 91 is 

K = + |6.458 x 10‘ 10 + 51.885 x 10" 6 00003 - cos (X + 156 . 0° )]| 1/2 (rad/sid. day) (92) 

(see Table 2, Case 2, and Figure 3). 


The Theory for 24-Hour Satellite Drift in a "J 33 -only” Perturbation Field 

Combining Equations 71, 2 and 67c causes Equation 71 to become for "j 33 -only" drift 


(*)’ 


2 A, 


= c 0 - 


— cos 3 (A-A 33 ) 


(93) 
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where 


[cos i s + l] 3 

A 33 = -540 7 7 2 ( R 0 / a s ) 3 J 33 g (rad/sid. day ) 2 


(94) 


The corresponding longitude acceleration during " j 33 -only drift" is given from Equation 66 as 

\ - A 33 sin 3(^ - A- 33 ) . (94a) 

Equation 94a establishes a one-dimensional potential energy for the longitude excursions of the 
24-hour satellite with respect to the J 33 perturbation. As before, this is given as the negative of 
the first integral of the right hand side of Equation 94a (ignoring the integration constant): 

A 

(P-*-fc 33 = “3 1 cos 3 (A. -^ 33 ) - (94b) 

(see Figure 9). 



Figure 9— 24-hour satellite longifude-driff potenfial 
energy and drift regimes due to J 33 . 


r . 2 A 33 

A = ± {( X o) 2+ ~3~ [cos 3(k 0 -\ 3 j) 


Since A 33 is always positive, the libratory 
longitude drift regime that can be established 
by j 33 for the 24-hour satellite has the form 
shown in Figure 9. Stable librations can take 
place only about the three longitudes A = a 33 ±60° 
and A = \ 33 + 180° , though momentarily no long 
term acceleration is experienced at A = \ and 
A = \ 33 ± 120 ° (see also Appendix A). 

If the initial longitude drift rate of the sat- 
ellite is A 0 with the ascending equator crossing 
at A 0 , then C 0 in Equation 93 is evaluated as in 
Equation 78 to yield for "J 33 -only" drift 

cos 3(A -A 31 )l\ 1/2 (rad/sid. day) . (95) 


(The sign of the root in Equation 95 is established from the sense of the drift with respect to the 
longitude where A = 0 for the libratory regime.) 

From Equations 72, 75 and 94, A 33 is 


2A 33 

— 3 — (i = 0°) = 122.9859 x 10 " 6 (rad/sid. day ) 2 


(96) 
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and 



(i 


60.0°) 


51.88468 x 10 6 (rad/sid. day) 2 


(97) 


For the case of the equatorial 24-hour satellite in Table 2, (Case 1) 


>^ 0 = -66.01067° 


and 


/; \ „ _ n 01066 X .99726961 2 _ . , . . . „ 

(^■o) “ I 57.29578 J - 3.442661 x 10 8 (rad/sid. day) 2 . (98) 

Thus substitution of Equations 72, 96 and 98 into Equation 95 gives the theoretical drift rate for 
the equatorial 24-hour satellite in Table 3, Case 1, as 

k = - |3.4427 x 10“ 8 + 122.99 x 10“ 6 [-. 00056 - cos 3(\ + 36 . 0° )]V /2 < rad / s id - da y ) • (99) 


See Table 3, Case 1 and Figure 4 for an evaluation of Equation 99 at the maximum longitude drift 
excursion of the satellite in the numerically computed trajectory. 

For the case of the 60.0° inclined- orbit, 24-hour satellite in Table 3, the theoretical drift 
formula corresponding to Equation 99 is 

k = " {7.9022 x 10~ 9 + 51. 885 x 10“ 6 [-. 00024 - cos 3(k + 36 . 0°)]| 1/2 ( rad/sid. day). (100) 

(see Table 3, Case 2 and Figure 4). 


The Theory for 24-Hour Satellite Drift in a "J 42 -only” Perturbation Field 

Combining Equations 71, 2 and 67d causes Equation 71 to become for " J 42 -only ,r drift 


(\> = C 0 -A 42 cos2(^-\ 42 ) 


where 


( 101 ) 


ISO/, 2 (R 0 /a,)« J 4 


( 1 + cos i ) 2 

4 


7 sin 2 i cos i(l + cos i) 

4 


( 102 ) 
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The corresponding longitude acceleration during "J 42 -only drift" is given from Equation 66 as 


A 4 2 s in 2 ^ 2) 


( 103 ) 


Equation 103 establishes a one-dimensional potential energy for the longitude excursions of the 
24-hour satellite with respect to the J 42 perturbation. As before, this energy function is given as 
the negative of the first integral of the right hand side of Equation 103 (ignoring the integration 
constant): 


(P - E ) ^42 ” COs2 ( k ~ X 42) 

(see Figures 10 and 11). 


(104) 



Figure 10—24-hour satellite longitude-drift potential 
energy and drift regimes due to J 42 (A 4 2 < O, -34. 1° 
< i < 34.1°, |i| > 80.25°). 


Figure 11—24-hour satellite longitude-drift potential 
energy and drift regimes due to J 42 (A 42 > O, 34.1° 
< lil < 80.25°). 


Depending on the inclination of the 24- hour satellite, two oppositely directed libratory longitude 
drift regimes are possible due to J 42 (also see "NONRESONANT INCLINATIONS FOR THE EARTH 
GRAVITY HARMONICS J 3l AND J 42 ). In the first, for |i| <34.1° and |i| >80.25° , stable long 
term librations can occur about \ = \ 42 . In the second regime, for 34. 1° < | i I <80. 25° , stable 
long-term librations can occur about A. = A 42 ±90° (see Figures 10 and 11 and Appendix A). At the 
critical inclinations i = ±34. 1° , ±80.25° and ±180.0° , there is no long term drift perturbation at 
any longitude due to J 42 , since A 42 = 0 for these inclinations. 

The Theory for 24-Hour Satellite Drift in a "J 44 -only” Perturbation Field 

Combining Equations 71, 2 and 67e causes Equation 71 to become for "j 44 -only" drift 




a 44 

C 0 ~ ~T COS 4 ( X " X 44) 


( 104 ) 
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where 


A 44 = -S040^(R 0 /a s )«J 22 5^ • (105) 

The corresponding longitude acceleration during "J 44 -only" drift is given from Equation 66 as 

k = A 44 sin 4(\ - \ 44 ) (106) 

Equation 106 establishes a one- dimensional potential energy for the longitude excursions of the 
24-hour satellite with respect to the J 44 perturbation. As before, this energy function is given as 
the negative of the first integral of the right hand side of Equation 106 (ignoring the integration 
constant): 


(P,E -^44 “ 4 cos4(^-^ 44 ) 

4 4 

(see Figure 12). 


(107) 


Since A 44 is always positive, the libratory 
longitude drift regime that can be established by 
J 44 for the 24-hour satellite has the form shown 
in Figure 12. Stable librations can take place 
only about the four longitudes k = \ 44 ±45° and 
k - \ 44 ±135°. Momentarily, no long term ac- 
celeration is experienced at the four longitudes 
K = j\ 44 , A = X 44 ±90° , and k = \ 44 + 180° (see 
Appendix A). 

It is noted that the actual longitude potential 
energy function of the earth for the 24- hour 
satellite will be a superposition of Figures 6 
through 12. The exact form of this total potential 
the actual longitude gravity constants of the earth 
1 and 10 that J 22 dominates the actual longitude d 
the earth. 



Figure 12—24-hour satellite longitude-drift potential 
energy and drift regimes due to J 44 . 


energy will depend on the relative proportions of 
It is evident from the studies in References 6, 
rift potential energy for the 24- hour satellite of 


Comparison of a Theoretical and Numerically Integrated 24-Hour Drift Trajectory in a 
Third Order Earth, Solar and Lunar Gravity Field 

In Table 4, Case 1 are tabulated the results of a numerically integrated 24- hour satellite 
trajectory in a third order earth, solar and lunar gravity field. The longitude earth field used 
is a recent weighted average estimate from W. M. Kaula (Appendix A). 
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The theoretical phase plane trajectory of the long term drift of a 24-hour satellite in a full 
earth field alone, is derivable from Equation 71. Let k Q be the initial small drift rate of such a 
satellite at a time when its mean longitude is at . Then from Equation 71 the constant C 0 is 
evaluated as 




cosm(K 0 -k nm ) F(i) n „ 


(108) 


Substitution of Equation 108 into Equation 71 gives the drift rate as a function of longitude: 


k = 


r f 1/2 

jw 2 + ^T 2Z F (i)nm [cosm(\-\ nm ) -cosm(\ 0 -\ nm )H 


(109) 


In terms of the A nm gravity-orbit constants, Equation 109 can be rewritten as 


k = 


~ 1/2 

|(\>) 2 + 7^ [cosm(*. 0 -X. nn ) -cosm(\-\ nm )]j 


( 110 ) 


The A nm f s for the trajectory in Table 4, Case 1 are the same as 71k. For these constants and the 
given initial conditions, the theoretical trajectory of Table 4, Case 1 is 


k 


^“4.2556 x 10 6 ^ y \ cos (rad/s id. day) 


dll) 


(see Table 4, Case 1). 

It is noted that the maximum deviation of the actual from the theoretical phase plane trajectory in 
this case is of the order of 0.5% and is due to long term effects of the sun and moon (see Figure 5a). 


"NONRESONANT” INCLINATIONS FOR THE EARTH GRAVITY HARMONICS J 3 i AND J 42 

Equations 67a-67e make it evident that at every inclination the total long term drift regime is 
due to different proportions (functionally independent) of the relevant gravity harmonics. In gen- 
eral, it can be noted that all the harmonics higher than the second have reduced influence with 
respect to the J 22 harmonic for increasing inclinations (at least for prograde orbits). For the 
"mixed tesserals" J 31 and J 42 , which are strongly latitude- sensitive, a large part of the smoothing- 
down effect is due to the latitude excursion of the inclined orbit. For the higher spatial frequency 
sectorial tesseral harmonics J 33 and J 44 , the effect of the daily longitude excursion is to smooth 
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down their orbit- averaged contributions with respect to that of the J 22 harmonic. (See the pictorial 
representations of the potential in Appendix A.) 

These facts make it possible, for example, to distinguish between the effects of the gravity 
harmonics j 22 and J 42 by utilizing two 24-hour satellites of reasonably different inclinations. 

With a single 24-hour satellite, this would not be possible because these harmonics have the same 
spatial frequency m. Of special interest for geodetic satellites are those inclinations for which 
the "tesserals," (j nm with n-meven, n t m) have zero inclination factors F(i) nm . Satellites 
having these special inclinations and periods commensurate with the rotation period of the earth, 
such as the 24-hour satellite, can have no secular perturbations of the semi- major axis due to that 
particular mixed tesseral harmonic. It is conjectured that there can always be found at least one 
"nonre sonant” inclination for each of these harmonics for which F(i) nm = 0. Secular perturbation 
is defined as a non- zero along track force which is orbit averaged over the synodic period of the 
satellite with respect to the earth’s rotation. Utilization of these special inclinations would sim- 
plify the difficult problem of separating the effects of the individual longitude harmonics in de- 
termining the exact shape of the geo id from observations on a limited number of earth satellites. 


"Nonresonant” Inclinations of J 31 

In Equation 67b let x = cos i. Then the condition for a nonresonant 24-hour satellite due to J 31 
is established when 

F(i) 31 = 0 = 1 2 ~ ~ ( 1 + 3x ) = 4 + 4x - 5 + 5x 2 + 15x 3 - 15x . (112) 

Equation 112 reduces to 

(x + 1) (l5x - 10X- l) = 0 . (113) 

The three roots of Equation 113 (giving the critical inclinations for J 31 ) are 

Xj = -1C i = ±180°) 

x 2 = .755(1 = ±41.0°) 

x 3 = -.089(1 = ±95.1°) * (114) 

"Nonresonant” Inclinations of i 42 

In Equation 67d let x = cos i . Then the condition for a nonresonant 24-hour satellite due to J 42 
is established when 

F( i ) 42 = 0 = (1+x) 2 - 7 (l - x 2 )(x) ( 1 + x) , (115) 
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Equation 115 reduces to 


(1+x) 2 (7x 2 -7x+l) - 0 (116) 

The four roots of Equation 116 (giving the critical inclinations for j 42 ) are 

Xj , x 2 = — 1( i = ±180°) , 

x 3 = . 828( i = ±34.1°) , 

x 4 = • 172( i = ±80.25°) . (117) 


DISCUSSION 

The exact theory of 24-hour satellite drift due to earth gravity presented in this report and 
summarized by Equations 66, 68 and 71 rests on three fundamental assumptions. The first is that 
the orbit of the satellite is always perfectly circular. It should be evident that any eccentricity in 
the orbit will give rise to greater errors in this theory for inclined satellites than for equatorial 
ones. For the equatorial satellite of up to moderate eccentricity (i.e., e < .2 ), the only change in- 
dicated in the simple theory would be to use the mean daily longitude position of the satellite as 
the characteristic longitude \ of the satellite. In fact, this is also the indicated theory improve- 
ment change for the equatorial satellite of appreciable drift rate (i.e. k > 2°/day) since the 
perturbation method relies fundamentally on orbit averaging the earth gravity effects with respect 
to geographical longitude and latitude. The situation with respect to eccentricity for inclined- 
orbit satellites is more severe. This is not merely because of the greater difficulty in locating 
the orbit averaged longitude location of the inclined- eccentric orbit 24-hour satellite. Difficult as 
this may be for perigee locations distant from the equator, the essential phenomenon which might 
cause distortion of the simple circular-inclined orbit theory applied to orbits with eccentricity is 
the possible introduction of significant energy perturbations to the orbit due to zonal earth gravity. 

W. M. Kaula* believes that long term effects on the semi-major axis of a satellite due to 
zonal gravity must be negligible to at least the second order in the eccentricity. Assuming con- 
tributions to orbit- averaged energy from zonal gravity proportional to the third power of the ec- 
centricity, it is easy to see that these will give rise to effects of the order of magnitude of the 
longitude effects only when the eccentricity is of the order of magnitude of 0.1. This is because 
the magnitude of the dominant zonal gravity force is of the order of 1000 times that of the strongest 
tesseral force at a given altitude. 

The second assumption is that the semi- major axis of the 24-hour satellite is unchanged in 
the drift from a perfectly synchronous orbit. In Reference 6 it was shown that over the widest 


♦Private communication. 
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possible longitude libration for a 24-hour satellite the expected maximum change of semi- major 
axis is only of the order of 0.1%. 

The third assumption (related to the second) is that the period of the orbit remains synchronous. 
This assumption enters the theory in two places. One is in the calculation of the orbit- averaged 
effect of the harmonics. These effects are assumed to be referable to a single constant longitude 
location over each orbit, an assumption valid only for an orbit with the earth synchronous period. 

The second place is in the assignment of the units ,T sidereal day" for the basic orbit-averaged drift 
Equations 9 and 11 and all subsequent drift equations. The proper time unit in these orbit averaged 
equations is the period of the orbit 2tt/co s , where o> s is the mean motion (mean angular rate) of the 
satellite in its orbit. LeW = \co s -co e | measure the difference of the actual mean motion of the 
satellite from the earth rate. Then it is evident that the error arising from the assumption of a 
constant earth period unit in the drift equations is of the order of w'/w . It can be shown that the 
first source of error in the calculation of the orbit averaged effects (referred to a single longitude) 
for the longitude harmonic J qp , in the case of the nonsynchronous satellite of mean motion a/ , is of 
order pa//w e . 

Thus to limit errors in the simple drift equations for the second order harmonic j 22 and the 
fourth J 42 to one part in a hundred, 2co'/u e < io -2 . Since the drift rate for the near- synchronous 
satellite is given approximately as 360 co'/co e degrees/day, the limiting drift rate with respect to 
these harmonics is 1.8 degrees/day. With respect to the third harmonic, j 33 , it is 1.2 degrees/day. 
With respect to the fourth order harmonic, J 44 , it is .9 degrees/day. With respect to the third 
order harmonic, J 31 , the limiting drift rate is 3.6 degrees/day. For drift rates in excess of these 
limits, the simple theory given in this report may still hold with reasonable accuracy if the longi- 
tude K is redefined as the mean geographic longitude traversed in the ground track over an orbit 
(i.e. the mean of successive ascending equator crossings for a nearly circular orbit). In addition, 
the time unit of the drift equations should be redefined in terms of o> s . 

Many numerical trajectories made in conjunction with geodetic analysis of the motions of the 
Syncom II satellite (References 6 and 7) have confirmed the adequacy, for these purposes, of the 
simple, "first-order" perturbation drift equations presented herein (at least for the J 22 effect) 
with unmodified \, for eccentricities up to .0012 and drift rates as high as 1 degree/day. 

CONCLUSIONS 

1. The first order long-term longitude drift of a 24- hour satellite due to earth gravity through 
fourth order has been shown to be strictly analogous to a mathematical pendulum for each of the 
relevant tesseral harmonics J 22 , J 31 , J 33 , J 42 and J 44 . 

2. Under the perturbation of the relevant gravity harmonic J nm , there are potentially m stable 
longitude positions around the equator about which long-term longitude librations of the 24-hour 
satellite can occur. 

3. The actual perturbed drift of the 24-hour earth satellite due to high-order earth gravity 
will, to first order, be a linear sum of the separate harmonic effects in a proportion depending on 
the magnitude of the actual gravity coefficients. 
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4. For the M real earth' 1 the long term drift of the 24- hour satellite is undoubtedly dominated 
by the second order tesseral effect (J 22 ). Terms of third order are probably next most influential 
on the 24- hour satellite, dominated by the effect of j 33 . 

5. The simple, first order pendulum drift theory due to earth gravity presented herein is 
valid for 24- hour circular orbit satellites of all inclinations and should give a good approximation 
to the actual long term drift for nearly synchronous period satellites having an eccentricity of up 
to order 0.1 and/or drift rates of up to the order of 1 degree/day. The theory can be expected to 
give the poorest accuracy for predicting drift in the immediate vicinity of the longitudes of the 
axes of the earth's elliptical equator under the additional perturbing influence of the sun and 
moon. 


6. The drift perturbation for each relevant gravity harmonic is strongest on the 24-hour 
equatorial satellite, and for the "tesserals," J 31 and J 42 , becomes zero at all mean longitude 
locations for certain critical orbit inclinations. 

7. For geodetic purposes, use of one or more inclined orbit 24-hour satellites (such as 
Syncom II) as well as equatorial 24-hour satellites gives promise of defining most efficiently 
gravity terms in the earth's potential to third and fourth order. 


(Manuscript received June 3, 1965) 
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Appendix A 


Earth’s Gravity Field as a Series of Spherical Harmonic Components. 

Comparison of Recent Longitude Gravity Fields 

The gravity potential used as the basis for this study is the exterior potential of the earth de- 
rived in Reference A-l for geocentric spherical coordinates referenced to the spin axis and center 
of mass of the earth. The infinite series of spherical harmonics is truncated after J 44 . The 
zonal harmonic constants J 20 , J 30 and j 40 used in the section on the comparison of the theory 
with numerically integrated trajectories are taken from Reference A- 2. 

The earth equatorial radius R 0 used in this study is 

R 0 = 6378. 165 km . 

The Gaussian gravity constant of the earth used is 

= 3.986032 x 10 s km 3 /sec 2 . 

The longitude harmonic constants used in the simulation studies are those of Kaula- Combined (1964) 
as listed in Table A-l. These correspond to the values shown on the "longitude geoids" below. The 
gravity potential of the earth (to fourth order, probably sufficient to account for all significant per- 
turbations on a 24-hour satellite) may be illustrated (following Reference A-3 with the zonal con- 
stants of Reference A-2) as follows, from Kaula- combined (1964): 
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The earth-gravity field (per unit test mass) whose potential is Equation A-l is given as the 
gradient of Equation A- 1 or 


„ dV E \ <9 V e 4> dV F 

F = fF r + KF X + <*F* = VV E = f TT + r + T1J * ( A '2) 

F r = 77 { _1 + ( R o/ r ) 2 [3/2J 20 (3 sin 2 (?i- 1) + 9J 22 cos 2 ?icos 2(>t- ^ 22 ) 

+ 2(R 0 /-) J 30 (5 sin 2 <t> ~ 3) (sin<£) + 6(R 0/ /r) J 3J (5sin 2 <£-l) cos <p cos (^~^ 31 ) 

+ 60(R 0/ /^) J 32 cos 2 (isin^cos 2(\- \ 32 ) +60(R 0 /r) J 33 cos 3 4> cos 3(\ - A- 33 ) 

+ 5/8(R 0/ /r) 2 J 40 (35 sin 4 0 - 30 sin 2 0 + 3) 

+ 25/2 (R o/r) 2 J 41 (7 sin 2 cb - 3) cos <f> sine fi cos (A. - \ 41 ) 

+ 75/ 2 (R 0 /r) 2 J 42 (7 sin 2 “ l) cos 2 $ cos 2 (^ - \ 42 ) 

+ 525(R 0/ /r) 2 J 43 cos 3 c£sin<£cos 3(\ - \ 43 ) + 525 (R 0 /r) 2 J 44 cos 4 c/> cos 4 (*"\ 4 )]} • (A- 3) 

r" 

F A = 77 ( R o/ r ) 2 |6J 22 .cos0sin2(\-A 22 ) + 3/2 (R 0 /r) J 31 [5sin 2 <£-l] sin(\-\ 3J ) 

+ 30(R„/r) J 32 cos $ sin 4> sin 2 (A. - X 32 ) + 45(R 0 /r) J 33 cos 2 sin 3 (\ - \ 33 ) 

+ 5/2(R 0 /-) 2 J 41 [7 sin 2 0-3] sin<£sin(\-\ 41 ) + 15(R Q /r) 2 J 42 (7sin 2 <£- l) cost^sin 2(\- \, 2 ) 

+ 315 (R 0 /r) 2 J 43 cos 2 sin sin 3 (a. - \ 43 ) 

+ 420(R 0 /r) 2 J 44 cos 3 <t> sin 4 (\ - \ 44 )| • (A-4) 

. v r 

F 0 “ ^7 ( R o/ r ) 2 V 3 J 2 o sin <£ cos <£ + 6J 22 cos <j> sin <£ cos 2 (\ - ^ 22 ) 

3/2(R 0 /r) J 30 (5sin 2 <£-l) cos 0 + 3/2 (R 0 / r ) J31 (l5 sin 2 cfc - ll) sin^cos 
+ 15(R 0/ /r) J 32 (3 sin 2 <£ - l) cos <f> cos 2 ~ ^ 32 ) 

+ 45(R 0 /r) J 33 cos 2 <£sinc£cos 3(\-\ 33 ) - 5/2 (R 0 /r) 2 J 40 (7 sin 2 <£- 3) sin<£cos<£ 

+ 5/2(R 0 /r) 2 J 41 (28 sin 4 <f>- 27 sin 2 0 + 3) cos (\-\ 41 ) 

+ 30(R 0 /r) 2 J 42 (7 sin 2 <j> - 4) cos sin <£> cos 2 (A- - \ 42 ) 

+ 105(R o /r^ 2 (4 sin 2 </> - l) cos 2 cos 3(X - ^ 43 ) 

+ 420(R 0/ /r) 2 J 44 cos 3 0sin0cos 4(\- \ 44 ) j* * (A-5) 
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The actual sea level surface of the earth is to be visualized through Equation A-l as a sphere of 
radius 6378 km, around which is superimposed the sum of the separate spherical harmonic de- 
viations illustrated. To these static gravity deviations, of course, must be added a centrifugal 
earth- rotation potential at the surface of the earth to get the true sea level surface (Refer- 
ence A-l). 

Table A-l gives recent estimates from a wide variety of sources of longitude gravity constants 
in the earth's potential. 
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Table A-l 


Longitude Coefficients in the Earth’s Gravity Potential 


{ V E = TF L L [ 1 "(V r ) 1 

k n = 2 m = 0 u 


P^Csin^) ] nm 



As Reported 1915-1964 2 . 


Longitude Geoid 7 

J 22 

*22 

J31 

*31 

J32 

*32 

J33 

*33 

J 4 1 

*4 1 

J42 

*42 

J43 

*43 

J44 

*44 

(1) Wagner (1964)3 

-1.7 x 10" 6 

-19* 















(2) Kaula- Combined (1964)^ 

-1.51 

-15.5 

-1.51 x 10" 6 

0.0° 

-.102 x 10“ 6 

0.0° 

- .149 x 10* 6 

22.8° 

- .465 x 10" 6 

-136.0° 

-.163 x 10" 6 

37.0° 

-.061 x 10" 6 

- 1.9° 

-.0053 x 10" 6 

35.8° 

(3) Izsak (1964)* 

-1.00 

-17.0 

- .934 

-15.5 

-.116 

19.0 

- .173 

38.0 

- .949 

-146.0 

-.074 

47.5 

-.024 

- 3.9 

-.0206 

25.3 

(4) Kaula (1964)* 

(5) Anderle and Oester- 

-1.77 

-18.2 

-2.12 

- 5.4 

-.379 

10.5 

- .105 

23.1 

- .263 

-239.0 

-.117 

42.3 

-.0473 

15.0 

-.0104 

14.5 

winter (1963)3 

-2.09 

-14.1 







- .773 

-140.6 

1 

M 

OO 

—3 

34.7 

-.163 

- 4.3 



(6) Guier (1963) 3 

-1.80 

-10.4 

-1.77 

6.3 

-.286 

- 2.6 

- .204 

24.1 

- .73 

-141.0 

-.273 

38.6 

-.0791 

- .7 

-.0102 

35,0 

(7) Kaula (Sept. 1963) 4 

-1.51 

-18.1 

-1.65 

5.3 

-.144 

46.4 

- .145 

15.8 

- .471 

-228.0 

-.078 

44.2 

-.0265 

22.6 

-.0038 

23.3 

(8) Izsak (July 1963)* 

-1.05 

-11.2 

-1.1 

3.2 

-.20 

-21.8 

- .14 

20.0 

- .43 

-132.1 

-.13 

37.0 

-.026 

11.5 

-.019 

14.8 

(9) Kaula (May 1963)* 

-1.4 

-21.5 

-1.6 

- 1.9 

-.15 

35.8 

- .156 

18.5 

- .53 

-233.7 

-.12 

44.5 

-.019 

10.7 

-.0038 

23.3 

(10) Cohen (May 1963) 3 

-2.08 

-14.1 







- .775 

-159.0 

-.288 

34.6 

-.162 

- 4.3 



(11) Kaula (Jan. 1963)* 

-1.62 

-21.4 

-1.81 

- 3.57 

-.145 

6.6 

- .112 

37.6 

- .479 

-245.5 

-.072 

47.7 

-.0988 

5.9 

-.0132 

28.4 

(12) Uotila (1962)5 

-1.52 

-36.5 

- .685 

-81.0 

-.409 

- 5.2 

- .398 

19.5 

- .238 

-127.0 

-.211 

14.6 

-.082 

- 9.3 

-.0142 

- 2.6 

(13) Kozai (Oct. 1962)* 

-1.2 

-26.4 

-1.9 

4.6 

-.14 

-16.8 

- .10 

42.6 

- .52 

-122.5 

-.062 

65.2 

-.035 

0.5 

-.031 

14.9 

(14) Newton (April 1962) 3 

-2.15 

-10.9 







-2.53 

-189.1 







(15) Newton (Jan. 1962)3 

-4.16 

-11.0 















(16) Kozai (June 1961)* 

-2.32 

-37.5 

-3.21 

22.0 

-.41 

31.0 

-1.91 

51.3 

- .262 

-196.5 

-.168 

54.0 

-.044 

-13.0 

-.054 

50.3 

(17) Kaula (June 1961) 6 

- .55 

-13.3 

-1.19 

20.6 

-.33 

- .9 

- .21 

22.6 

- .617 

-166.0 

-.14 

21.1 

-.031 

- .5 

-.008 

26.4 

(18) Izsak (Jan. 1961)* 

-5.35 

-33.2 















(19) Kaula (1961)5 

-1.68 

-38.5 







-1.15 

- 13.0 







(20) Krassowski (1961?) 

-5.53 

15.0 















(21) Kaula (1959) 5 

- .62 

-20.9 

- .98 

55.4 

-.11 

13.3 

- .19 

14.3 

- .46 

-132.3 

-.081 

48.6 

-.01 

-30.0 

-.02 

22.5 

(22) Jeffreys (1959) 5 

-4.17 

0.0 















(23) Uotila (1957) 5 

-3.5 

- 6.0 















(24) Zhongolovitch (1957) 5 

-5.95 

- 7.7 

-2.21 

-25.7 

-.628 

-26.4 

- .54 

13.0 

- .78 

-149.1 

-.080 

45.0 

-.051 

- 3.8 

-.0224 

15.9 

(25) Subbotin (1949) 5 

-5.5 
















(26) Niskanen (1945) 5 

-7.67 

- 4.0 















(27) Jeffreys (1942) 5 

-4.1 

0.0 

-2.1 

0.0 

-.66 

0.0 

- .24 

33.3 









(28) Heiskanen (1928) 5 

-6.34 

0.0 















(29) Heiskanen (1924)5 

-9.0 

18.0 















(30) Helmert (1915)5 

-6.0 

-17.0 
















r is the radial distance of the field point to the center of mass of the earth, n the earth’s Gaussian gravity constant * 3.9860 x 10 20 cmVsec 2 , R Q the mean equatorial radius of the earth = 6378.2 km. $ is the geocentric latitude of the field 

point. A is the geographic longitude of the field point. J 21 = 0, since the polar axis is very nearly a principal axis of inertia for the earth. p“ (sin <f>) ~ cos'" 4> V T_. sin""" 9 ' 21 <t>, where K is the integer part of (n -m)/2 and 
T (-I) t (2n ~ 2t)l L 


nmt _ (See Kaula, 1964). The longitude coefficients are those for which m /o. 

2 tf (n t)* (o 01 2t)l 

The Jnm’* and A nm’ s in *** table > e * ce P r in one or two instances, have been converted from the original author’s set of gravity coefficients. The blanks indicate the author did nor consider that particular harmonic in fitting an earth potential 


to the observed data. In one or two instances, noted below, the author reported tesseral coefficients to higher order than the fourth. 
^Satellite - Doppler geoid 
* Satellite-camera geoid. 

^ Surface •gravimetric geoid. 


^Combined astro -geodetic, gravimetric and satellite geoid. 

7 Detailed information on references for the geoids listed below given on following page. 



References for Geoids Listed in Table A-l 


Geoid References 

(1) "Determination of the Elliptic ity of the Earth’s Equator from Observations on the Drift 

of the Syncom II Satellite” TN D-2759, May 1965. 

Private communication from W. M. Kaula (July 1964). 

Private communication from W, M. Kaula (July 1964). 

Private communication from W. M. Kaula (July 1964). 

Guier, W. H., and Newton, R. R., "Non- Zonal Harmonic Coefficients of the Geopotential 
from Satellite Doppler Data," Silver Spring, Maryland: Johns Hopkins University 
Applied Physics Laboratory, TG-520, November 1963. 

(7) Kaula, W. M., "Improved Geodetic Results from Camera Observations of Satellites," 

J. Geophys, Res. 68(18):5183-5190, September 15, 1963. 

Izsak, I. G., "Tesseral Harmonics in the Geopotential," Nature , 199(4889): 137-139, 

July 13, 1963. 

Private communication from W. M. Kaula (May 1963). 

Kaula, W. M., "Tesseral Harmonics of the Gravitational Field and Geodetic Datum 
Shifts Derived from Camera Observations of Satellites," J. Geophys. Re s. 68(2):473-484, 
January 15, 1963. 

Private Communication from W. M. Kaula (July 1964). 

Private communication from Y. Kozai (October 1962). 

Private communication from Newton (April 1962). 

Newton, R. R., "Ellipticity of the Equator Deduced from the Motion of Transit 4A," 

J. Geophys. Res. 67(1):415-416, January 1962. 

Kozai, Y., "Tesseral Harmonics of the Gravitational Potential of the Earth as Derived 
from Satellite Motions," Astron. J. 66(7):355-358, September 1961. 

(17) Kaula, W. M., "A Geoid and World Geodetic System Based on a Combination of 

Gravimetric, Astrogeodetic and Satellite Data," J. Geophys. Res. 66(6): 1799-1811, 

June 1961. 

Izsak, I. G., "A Determination of the Ellipticity of the Earth’s Equator from the Motion 
of Two Satellites," Astron. J. 66(5):226-229, June 1961. 

Kaula, W. M., "Analysis of Satellite Observations for Longitudinal Variations of the 
Gravitational Field," In: Space Resea rch II, (H. C. van de Hulst, C. de Jager, and 
A. F. Moore, eds.): 360-372, Amsterdam: North-Holland Publishing Co., 1961. 

(20) "Passive Dynamics in the Space Flight," Bureau of Naval Weapons Paper by J. D. 
Nicolaides and M. M. Macomber, 1962. 

(21) Kaula, W. M., "Statistical and Harmonic Analysis of Gravity," J. Geophys. Res. 64(12): 
2401-2421, December 1959. 

(22) Jeffreys, H., "The Earth: its Origin, History, and Physical Constitution, 4th ed., 

London: Cambridge University Press, 1959:127-194. 

(23) , (26), (28), (29), (30) Heiskanen, W. A., and Vening-Meinesz, F. A., The Earth and it s Gravity Field , New 

York: McGraw'-Hill Book Co., 1958:79 

(24) Zhongolovich, I. D., "Potential of the Terrestrial Attraction," Akad, Nauk SSSR Inst. 
Theoret. Astron. 6(8):505-523, 1957 (In Russian). 

Subbotin, M. F., Kurs Nebesnoi M ekhaniki, Vol. 3, 2d ed., Moscow: Gostekhizdat, 1949. 
In: Monthly Notices R oyal Astronom. Soc. Geophys. S uppl,, 5(55): 1942. 


(18) 

(19) 


(12) 

(13) 

(14) 

(15) 

(16) 


( 8 ) 

( 9 ) 

( 11 ) 


(2) 

(3), (4) 

( 5 ) 

(6) , (10), (14) 


(25) 

(27) 
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Appendix B 


Symbols Used 

A nm The long term 24-hour satellite longitude- drift driving function of the gravity harmonic J nm 

a, a s The instantaneous semi- major axis, and a near- synchronous semi- major axis of the 24- 
hour earth satellite assumed to be constant in a particular drift regime 

F A gravity force per unit mass, acting on a 24- hour satellite 

F t The tangential (or circumferential) component of a gravity force acting on a 24- hour near 
circular orbit satellite 

F An orbit-averaged gravity force 

F(i) nin The inclination factor of the longitude drift driving function A nm for the 24-hour satellite 

g s The radial acceleration of the earth's gravity field at the altitude of the synchronous 
satellite (^.7355 ft./sec. 2 ) 

i, i s The inclination of the orbit plane with respect to the equator and the inclination of the 
near synchronous satellite's orbit plane assumed to be constant in a particular drift 
regime 

Jnm’ Spherical harmonic constants (order n, power m) expressing magnitude and phase angle 

of that component of the earth's gravity potential 

R 0 The mean equatorial radius of the earth (~6378.2 km) 

a The angle between a meridian plane through the satellite's position and the orbit plane 
( ) 0 Time zero quantity of the argument ( ). 

6 The argument from the ascending node to the 24-hour satellite’s position (defined in Figure 1) 

V, r, ^ Geographic longitude, geocentric radius and geocentric latitude of the 24-hour satellite 
position. In the text ^ also refers to the daily longitude position of the ascending 
equator crossing, or the orbit- averaged, mean daily longitude position of the 24- hour 
satellite 

m e The earth's Gaussian gravity constant (3.98603 x 10 5 km 3 /sec 2 ) 
a;, cl> c A circular frequency, and the earth’s sidereal rotation rate (.7292115 x 10~ 4 rad. /sec.) 
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